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Abstract 

In this paper we review a part of the approaches that have been considered to explain 
the extraordinary discovery of the late time acceleration of the Universe. We discuss 
the arguments that have led physicists and astronomers to accept dark energy as the 
current preferable candidate to explain the acceleration. We highlight the problems 
and the attempts to overcome the difficulties related to such a component. We also 
consider alternative theories capable of explaining the acceleration of the Universe, 
such as modification of gravity. We compare the two approaches and point out the 
observational consequences, reaching the sad but foresightful conclusion that we will 
not be able to distinguish between a Universe filled by dark energy or a Universe where 
gravity is different from General Relativity. We review the present observations and 
discuss the future experiments that will help us to learn more about our Universe. This 
is not intended to be a complete list of all the dark energy models but this paper should 
be seen as a review on the phenomena responsible for the acceleration. 

Moreover, in a landscape of hardly compelling theories, it is an important task to 
build simple measurable parameters useful for future experiments that will help us to 
understand more about the evolution of the Universe. 
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1 Introduction 

Modern cosmology was born in the 1920s as one of the first applications of General Relativity 
for studying the properties of the cosmos on large scales. If the traditional astrophysics and 
astronomy were studying single objects (e.g. stars) or clusters of such objects (e.g. galaxies), 
characterizing the local properties of the Universe, with a more sophisticated instrumentation 
it was possible to observe deeper in the sky and at larger scales. More generally, if the global 
properties of the cosmos are to be described, we refer to the Universe as the first object 
under investigation by cosmology. 
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Modern cosmology has a precise birth-date: the discovery of Cepheids and ordinary stars 
in Nebulae by E. P. Hubble. After he recognized that most of the nebulae were galaxies, 
Hubble also pointed out that they were moving away from each other. Such expansion, 
currently named the Hubble flow, has been confirmed over the years by the whole set of 
data we now possess. Even if Hubble's intuition was extraordinary, his data were too noisy. 
At the distances where he pretended to see an average expansion, the motion of galaxies is 
still dominated by their peculiar motions. Therefore local measurements of moving objects 
are not accurate to prove the total expansion of the Universe. It is at scales larger than 
the galactic scale in which the dynamic evolution is mostly due to pure gravity. Galaxies, 
therefore, are the occupants of the Universe whose rules are set by general relativity and 
their average distances are increasing with time. 

An expanding Universe filled only by ordinary matter (e.g. matter whose stars are formed) 
will at some point slow down due to gravitational attraction. In 1998, astronomers, observing 
supernovae, realized that their apparent luminosities were reduced as an effect of more red- 
shifting; astronomers pointed out the Universe is not only in phase of expansion but this 
expansion is accelerating. Although the supernova observations, see Refs. |T] and [5] has led 
to the general acceptance that the expansion of the Universe is accelerating, we are still far 
from a consensus on what is responsible for the acceleration. Although there is a name for 
the phenomenon, dark energy, there is as of yet no convincing physical explanation. 

Nonetheless, other recent observations have put serious doubts on the completeness and 
understanding of the cosmological models, the theoretical and experimental considerations 
are: 

• Primordial nucleosynthesis: predicts the baryonic abundance which is in agreement 
with the observations but it is not enough to produce the observed gravitational effects. 
This has led people to believe that a new form of matter, which is not (or weakly) 
coupled to the electro-magnetic field, has to be added in the Universe, known as Dark 
Matter. 

• Inflation: (needed to solve the curvature and horizon problems) requires an extension 
of standard model of particle physics. 

• Supernovae type la (SNIa): show the Universe is in a phase of accelerated expansion. 

• Large Scale Structure (LSS): provides another test for the existence of dark energy. It 
also shows there is a lack of matter. Therefore, this is also another evidence of dark 
matter. 

• Cosmic Microwave Background Anisotropics (CMB): suggest the total energy density 
in the Universe has to be equal to a particular critical value. 

The result of these observations has built a scenario in which the Universe is in a phase 
of accelerated expansion (SNIa) and its matter density (LSS) is not sufficient to set the 
curvature equal to zero (CMB). 

This review is organized as follow: we first introduce the basic equations for studying the 
dynamic of the Universe, section [5] and |31 in section U] we report in more detail the Hubble's 
flow and definition of cosmic distances. 

The observational evidence that has led cosmologists to accept dark energy are reported 
in section [SJ in section [5] we list the most studied models for dark energy whereas in section 
[7| are shown some of the models of modified gravity that are able to lead the Universe to a 
phase of accelerated expansion. In section |5] we compare the two different approaches, i.e. 
dark energy (or equivalently to add a new form of matter in the Universe) and modified 
gravity (or revising the laws of gravity). 
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In section [5] we present a list of possible test that a model has to pass in order to be 
considered a good dark energy candidate. 

In section [TU] we relate theory to observations and we show how phenomenological pa- 
rameters can be helpful to distinguish, say, a cosmological constant from a dynamical dark 
energy or the latter from some form of modified gravity. We end with the new observational 
tests in section [TTl 



2 Basic equations 

General Relativity theory has led to an increasing development in the understanding of our 
Universe. The most interesting aspect of this theory is that the dynamics of every object 
can be expressed by a single tensorial equation: 

_8ttG 

known as Einstein equation. Here G^^ represents Einstein tensor, which defines the geometry 
of the Universe, G is Newton's constant, T^,y denotes the energy-momentum tensor, which 
describes the matter and the energy content of the Universe and c is the speed of light. In 
other words, Einstein equation relates the geometry of the Universe to its content. 

Einstein equations are in general complicated non-linear equations and it is impossible to 
solve them analytically and only few exact, analytical solutions are known, involving strong 
symmetry constraints; however under the conditions of homogeneity and isotropy, which is 
approximately true at large scales, the most general line of element for an expanding Universe 
is the Friedmann-Lemaitre- Robertson- Walker metric (FLRW): 



ds2 = _di2 +a2(^) 



1 - Kr^ 



(2) 



where a (t) is the scale factor and represents the evolution in time of all the physical scales 
(every scale £o will change in time as £ (t) = a (t) £q); K is the curvature parameter (it can 
be —1, 1 or meaning an open, closed or fiat Universe respectively); t is the cosmic time 
and the coordinates r, and (j) are spherical comoving coordinates. 

In the FLRW spacetime the energy-momentum tensor takes the simple perfect fluid form: 

T'"' = (p+p)u'V+pg^'' (3) 

where u'^ is the four velocity of the fluid and p and p are the energy density and the pressure 
of a fluid, respectively. 

Einstein equation gives rise then to the Friedmann equations: 

. K 8ttG 

4vrG , „ , 
= ^(P + 3p) (5) 



where the dot denotes the derivative with respect to the cosmic time t and H = a/a is called 
the Hubble parameter which describes the expansion rate of the Universe. The energy density 
p should be thought as a sum running over all the forms of matter that fill the Universe, the 
same discussion has to be made for the term {p + 3p). 

In both the Friedmann equations the left hand side is given by the left hand side of 
Einstein equation, in other terms the geometry of a homogeneous and isotropic fluid is fully 
characterized by the scale factor a {t) , whereas the right hand side of Friedmann equations 
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come directly from the energy momentum tensor T^'^ . Eq. ([3]). hence geometry is related to 
matter. 

The Einstein tensor satisfies Bianchi's identity: 

G^,, = (6) 
where represents the covariant derivative; from Eq. ([1]) follows automatically: 

T|^■,^. = (7) 
which gives the continuity equation of a fluid in the FLRW background Universe: 

p + iH{p + p)=Q. (8) 

Let us now rewrite Eq. ^ as: 

^W"l = -|^ (9) 
[aH] 

where D.{t) = p (t) / pc (t) is the dimensionless density parameter and pc (t) — '3H^ /SttG 
is the critical energy density parameter. Then, the matter content determines the spatial 
geometry of the Universe, i.e.: 

V, > 1 -> K = +1 ^ closed Universe (10) 
n = l is: = ^ fiat Universe (11) 

fl < 1 K = -I ^ open Universe. (12) 

Recent measurements of Cosmic Microwave Background Radiation (CMBR) have shown that 
the Universe is spatially flat (i7 = 1 or iiT = 0) to within a few percent, see Ref. [3] . This is 
also a result that comes naturally from inflation in the early Universe, see Ref. [18j. So we 
assume hereinafter, if not else specified, that K is equal to 0. 

3 The evolution of the Universe with a barotropic perfect 
fluid 

Let us now consider the case in which the Universe is filled only by one single perfect 
barotropicl^ fluid with an equation of state parameter: 

w^-. (13) 
P 

If w is constant one can analytically find the evolution for p and the scale factor a (t) by 
solving the Eqs. Q and (|5]): 

p oc a-3(i+«') (14) 
a(t) cx (t-to)^HTT^ (15) 

where to is a constant. An analytical expression for the Hubble parameter can be found: 

H{t)^ ^ (16) 

The above solutions, Eqs. and are valid for all the value of uu except for w = — 1. 
When the parameter of equation of state w is —1, i.e. p — ~~p it follows from Eq. (fT5)) that 



We define a fluid to be barotropic if the pressure p depends strictly only on the energy density p: p = p (p) 
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the energy density p is constant; this corresponds to the cosmological constant, denoted by 
A. In this case the Hubble parameter H is also constant and the scale factor evolves as: 



a(t)=exp(m) (17) 

which corresponds to the de-Sitter Universe. It is clear from Eq. (11711 that a cosmological 
constant cannot be an appropriate candidate for inflation in the early Universe because 
otherwise inflation would not have ended. In turn a cosmological constant can be a good 
candidate for the dark energy today; for the moment we assume the cosmological constant 
to be the only model for dark energy. 

The parameter of equation of state w fully characterizes the fluid at background level, 
for instance: 

Rad. : w — ^-, a {t) oc {t — to)^^'^ , p cx (18) 

Dust : w = 0, a{t)(x{t-'tof^^ , p(xa^^ (19) 
A : w ~ a (t) (X exp (Ht) , p ~ const. (20) 

Both cases, radiation and dust, correspond to a decelerated Universe; this can be easily seen 
if we insert the corresponding value of the parameter of equation of state in the Eq. ([5]), 
which gives d < 0. In order to have cosmic acceleration, we need d to be positive, in other 
words we require: 

w < -1/3, p < -p/3. (21) 

This is clearly an upper bound for the parameter of equation of state; we can say that: the 
more the w is negative the bigger d is. 

In order to explain the late phase acceleration in the Universe we require an exotic form 
of energy with negative pressur^l satisfying Eq. (|2ip . 

This concept is completely new and it comes directly from Einstein equations because of 
relativistic effects (absent in Newtonian gravity). To show this we can consider a homoge- 
neous sphere with radius r and energy density p in which a point particle with mass m is 
moving. In Newtonian gravity, the equation of motion for such a particle is given by: 



mr = — - 



leading to: 

- = 23 

r 3 

which looks like Eq. ([5]) except for the missing pressure term. Then, we need a large negative 
pressure, to counter-balance the gravitational attraction, in order to accelerate the Universe. 
As shown before, Newtonian gravity can only decelerate the expansion of the Universe. 



4 Hubble's flow and Cosmic Distances 
4.1 Hubble's flow 

As mentioned before, in 1929 E. P. Hubble, see Ref. [19], observed that the spectral lines of 
the chemical elements in the galaxies were shifted from their normal position, they manifested 
a shift towards the red part of the electro- magnetic spectrum (hence the name redshift), more 
precisely towards lower frequencies with respect to those obtained in the laboratories. Such 

•^The energy density is assumed to be positive. 
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a characteristic, analogous to the Doppler effeclQ. is expressed in terms of parameter z given 
by: 

z = ^ (24) 

where Aq and Ae are the observed and emitted wavelengths, respectively. He concluded that 
objects were moving away from us and he also derived, for small redshift, a direct relation 
between the redshift and the distance of an object: 

z = HqD. (25) 

The redshift increases with distance or equivalently as we go back to the past the redshift 
increases. 

Here Hq is the present value of the Hubble parameter called Hubble constant, which is 
usually written as: 

Ho = 100 h Km sec-^Mpc"^ (26) 

where h describes the uncertainty of the Hubble parameter; measurements constraint this 
value to be, see Ref. |20] : 

/i = 0.72 ± 0.08. (27) 

Moreover, an object at the distance D will show a recession velocity proportional to its 
distance; to show this, let us consider two objects separated by a comoving distance ro, than 
their physical distance will be R{t) = a (t) ro , if we derive the latter we find: 

Vr^dro^-R (t) = HE (t) (28) 
a 

thus the relative velocity of the two objects is independent on the position of the observer. 
4.2 Cosmic distances 

An important concept related to observational tools in an expanding background is associated 
with the definition of distance. In fact most of the evidence in dark energy comes from 
distance measurements. Therefore, in order to discuss the observational constraints, we need 
to introduce the cosmic distances in the FLRW Universe. Let us rewrite the metric ^ in 
the following form: 

ds2 = -df + a {t) \dx^ + (Sk {x) f (29) 



where dfl — d9^ + sin 9'^dd)^ and where: 



{sin(x) K = +l 
X K = (30) 

sinh(x) K = -l 

4.2.1 Comoving Distance 

Let us consider the metric given by Eq. ([29l) . we can evaluate the physical distance along a 
null geodesic simply setting: ds^ = —c^dt^ + c? (t) dx^ = 0; that all light rays travel along 
null geodesic. Let us suppose that light is emitted by a source at redshift z (corresponding to 
Xe)) at time and reaches the observer situated at z = (corresponding to x = 0) at time 
to\ the comoving distance can be evaluated simply by integrating the null geodesic equation: 

= - r -^di = — (31) 



*But it is not. Firstly because, at very large distances we would have had objects moving faster than the 
speed of light. Secondly, in General Relativity the presence of matter induces a gravitational redshift. 
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4.2.2 Luminosity Distance 

Another way of defining a distance is ttirougli tlie luminosity of tlie objects. The luminosity 
distance plays an important role in observational cosmology. Let us consider a source with 
an intrinsic luminosity Lg at some distance ds , the flux we observe will be which is defined 
as amount of energy per unit surface, that is = £o/j the area of a sphere at z = in a 
curved space time is — 47r {aoSx (x)) ■ The luminosity distance can be written as: 

Dl^iaoSK{x)f^- (32) 

'-0 

The luminosity of an object is defined as the amount of energy emitted per unit time. If 
we denote by AEe the energy of light emitted by an object in the time interval Ate with 
absolute luminosity £s and by AEq the energy that reaches the observer, we have: 

C. - H (33, 

Since the energy of the photon is inversely proportional to its wavelength A we have that: 

^^^-<>-) 

where we have used Eq. ([M)) . 

Similar considerations can be made for the time interval, so we have that Ate/Ato — 
Ae/Ao = (1 + z). Hence we find: 

AEeAto 2 

Liserting Eq. ([55)1 into Eq. (|32p the luminosity distance reads: 

= ao5if + (37) 

4.2.3 Angular Diameter Distance 

The angular diameter distance is the distance under which an object of actual size Ace is 
seen under an angle A6, that is: 

D.^^. ^ ^ ^ (38) 

The size Ax of a source lying on the surface of a sphere or radius x at the time ti with the 
observer at the center is given, in FLRW metric, by: 

Ax ^ a{ti)SKix)^S- (39) 
The angular diameter distance will become: 

DA = a{t,)SK{x}^T^SK{x) (40) 
1 + z 

where we have used the relation z + 1 ~ oo/a (ii). Comparing the last one with Eq. p7p we 
find the relation, called distance duality: 

DL = {l + zfDA. (41) 

The distance duality relation holds for any metric (not just FLRW) in which photons travel 
on null geodesic and it depends crucially on photon conservatiorH, see Refs. [211 122). 

^For instance if photons are scattered from dust or free electrons, the luminosity distance will be 
reduced but the angular diameter distance is basically unaffected, as a consequence the relation H41| l does no 
longer hold. 
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5 Observational evidence for dark energy 



The existence of dark energy is supported by many observations. It was in 1998 when two 
groups, see Refs. \X\ and [2J, firstly and independently reported that the Universe is a phase 
of accelerated expansion by observing the luminosity of Supernovae Type la (SN la) which 
are huge explosions of some stars. This event is extremely luminous and can be used as a 
test in observational cosmology. The Supernovae are classified according to their spectral 
characteristics. The most useful in cosmology are those called Type la because they are poor 
in hydro and the mechanism generating this explosion is assumed to be known. The 
explosion of Type la occurs when the mass of a white dwarf in a binary system (usually 
its companion is a red giant star) exceeds the Chandrasekhar limit by absorbing the gas 
from the other star. The general belief is that SN la are formed in the same way in all the 
Universe, which means they have a common absolute luminosity. Hence SN la are considered 
to be standard candles by which the luminosity distance can be determined. In particular, 
the light curve of a SN la is correlated with its peak luminosity to a precision of about 7%, 
in other words, the brighter the SN the longer it takes to vanish, see Ref. [23;. 

Other observations, apart from SN la, support the evidence for dark energy and these 

are: 

• Age of the Universe compared to the oldest star. 

• Cosmic Microwave Background (CMB). 

• Baryon Acoustic Oscillation (BAO). 

• Large Scale Structure (LSS). 

5.1 Constraints from SN la 

In astronomy the magnitude is often used as quantity to categorize sources; also in this case 
we have the absolute magnitude A4 which is related to the logarithm of luminosity and the 
apparent magnitude m which is related to the flux. Absolute and apparent magnitude are 
related to the luminosity distance via the relation: 



Since the absolute magnitude A4 is known a priori (under the assumption of standard candle) 
then the luminosity distance can be derived from Eq. (|42l) once the relative magnitude m is 
measured. The redshift of the corresponding source can be obtained simply by measuring 
the frequency of the photons. 

Let us consider now, as an example, a FLRW flat background {Sk (x) — x)j the lumi- 
nosity distance Eq. (|37| reads: 



where E (z) ~ H (z) /Hq. The Hubble parameter can be expressed in terms of luminosity 
distance: 



Once the luminosity distance is measured by observations then we can determine the expan- 
sion rate H (z) and this is crucial to understand our Universe because we know the Hubble 

^It is the most abundant element in the Universe. 




(42) 




(43) 




(44) 



9 



parameter depends on all the components present in the Universe. The energy density in 
Eq. (HI has to be thought as the sum over all the components, namely relativistic matter, 
non-relativistic matter, cosmological constant (and eventually other stuff): 

p-j:po. (^)"''^^' =j:poAi+^r'^'"'- (45) 

Using Eq. the Hubble parameter reads: 

H^z)^H^Y.^oA^ + ^^f""^ (46) 

i 

where Oo.i = Po,i/po.c = '^H^po^i/SnG, where we have used the critical energy density defined 
previously. 

As an example, let us consider a Universe filled only by non-relativistic matter w — and 
cosmological constant w = —1 (relativistic matter plays a small role in the late time and can 
be neglected), the Hubble parameter and the luminosity distance become: 

(z) = Hi b„,o (1 + zf + l^A^ol (47) 



c 







dz' 



Dl{z) = 7r(i + ^)/ (48) 



In Fig. ([T]) we plot the luminosity distance Eq. (|48| . We remind the reader that we are 
working in a flat FLRW background which implies that fim.o + ^^A.o = 1- As we can see 
the luminosity distance increases with the value of the cosmological constant. Fig. ^ il- 
lustrates the observational values of the luminosity distance versus redshift z together with 
the theoretical curves derived from Eq. pS)) . The plot clearly shows that a Universe with 
only matter — 1) does not fit to the data. Perlmutter et al. showed that cosmolog- 

ical constant is present at the 99% confidence level, giving the value of the matter energy 
density Vtmo = 0.28lQQg (Icr statistical) in the flat Universe with a cosmological constant 
which has to be roughly the 70% of the total energy density of the Universe. Since at low 
redshift the models are indistinguishable and moreover observational data suffers from sta- 
tistical and systematic errors, one cannot conclude that the Universe is accelerating from a 
small data set. Recent experiments as the Super Nova Legacy Survey (SNLS), Hubble Space 
Telescope (HST) and Equation of State: SupErNovae trace Cosmic Expansion (ESSENCE), 
see Refs. HSl HBj respectively, have measured the luminosity distances of SN la at higher 
redshift giving tighter constraints on the dark energy parameters and conflrming the trend 
that the Universe is in a phase of accelerated expansion. 

5.2 Age of the Universe 

The inverse of the Hubble parameter gives roughly the age of the Universe. Once an ex- 
pression for the Hubble parameter is taken then one can find the corresponding age of the 
Universe and compare it with the oldest stars which gives a lower age boundary. Let us con- 
sider the Hubble parameter given in Eq. (|46ll in a FLRW Universe and taking into account 
all the component, i.e. relativistic and non-relativistic matter and dark energy: 

H (z) = i7o [^ro (1 + zf + (1 + zf + noE,, (1 + ^)3(i+---)] (49) 

where we have not made any assumption on the nature of the dark energy yet. Using the 
relation dt = — dz/ [(1 + z) H] we can find the age of the Universe to be: 

r°° dz 
Jo (l + z)i?(z) 
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Figure 1: The luminosity distance Dl for a flat cosmology with two components, matter 
Wm = and a cosmological constant wa = — 1 for different values of • 




Figure 2: The effective apparent magnitude ms for 42 high-redshift SN la from the Supernova 
Cosmology Project, and 18 low-redshift SN la from the Calan/Tololo Supernova Survey, after 
correcting both sets for the SN la lightcurve width-luminosity relation. The inner error bars 
show the uncertainty due to measurement errors, while the outer error bars show the total 
uncertainty when the intrinsic luminosity dispersion, 0.17 mag, of lightcurve- width-corrected 
Type la supernovae is added in quadrature. The horizontal error bars represent the assigned 
peculiar velocity uncertainty of 300 km s~^. The solid curves are the theoretical to^°'^*'™(z) 
for a range of cosmological models with zero cosmological constant: (f2M,f^A) — (0,0) on 
top, (1,0) in the middle and (2,0) on bottom. The dashed curves are for a range of flat 
cosmological models: (Qmi^^a) = (0,1) on top, (0.5,0.5) second from top, (1,0) third from 
the top, and (1.5,-0.5) on the bottom. From Ref. 2. 
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The integral in Eq. (|50p is too complicated to solve analytically and we are forced to make 
some assumptions. The contribution of the radiation term in Eq. (|50p can be neglected since 
the radiation dominated era is much shorter than the total age of the Universe, or to put it 
differently, the integral in Eq. ([50)1 is not hardly affected for redshift z > 1000. Moreover, 
^Iro is of the order of 10"'*, then it is a good approximation to neglect radiation. Let us now 
consider the case of a cosmological constant [w = —1), then the age of the Universe is: 



1 



dz 



(1 + z) n^^ {1 + zf + noEo 



1/2 ■ 



(51) 



integrating the last one, we have: 



to 



1 



3Ho yr 



nr. 



:log 



(52) 



where we used flmo + ^dEq — 1- In the limit Udeq ^ we find: 



to 



3Ho 



(53) 



A Universe without cosmological constant would have an age within the range 8.2 Gyr < to < 
10.2 Gyr, where we have used h = 0.72. However, different groups have estimated the age of 
the globular cluster, see Refs. |i27 l [28 l [29 ] . to be within the range of 12.7 Gyr < to < 13.5 Gyr. 
In most cases the globular clusters seem to be older than 11 Gyr (upper bound given by 
Eq. ([55]) ). being inconsistent with the age of a matter dominated Universe. The age of the 
Universe can be increased by taking into account the dark energy with an equation of state 
close to —1 as shown in Fig. ([3]). 

Also in an open Universe flxg > it is possible to increase the cosmic age but it will 
be still insufficient to bring the cosmic age up to 13Gyr. As an example, let us consider a 
Universe with + ^Kq = 1 ^^d zero cosmological constant, then the age of the Universe 
will be: 



to 



Hn 1 — 0„ 



1 



n 



mo 



loa 



1- VT^o; 



1 + v/r^a; 



(54) 



when ilmo — 1 recover Eq. (|53p ; in the limit of ^Iko — 1 then the age of the Universe is to — 
1/Ho — 13Gyr. However the curvature Hko has been constrained by WMAP measurements, 
see Ref. [3J, to be much smaller than unity. An open Universe without dark energy is 
inconsistent with the age of the oldest stars present in the Universe. The flat Universe with 
the cosmological constant is consistent with WMAP constraints for a value of matter density 
0.271 < Qmo < 0-289, see Fig. ©. 



5.3 Cosmic microwave background 

During the last decade an avalanche of experiments, together with the most recent WMAP 
satellite, have measured the small and large angular temperature fluctuations of the Cosmic 
Microwave Background Radiation. Such measurements have detected a series of acoustic 
peaks in the anisotropy power spectrum conflrming the earlier predictions about the evolution 
of pressure waves in the primordial photon-baryon plasma. The specific features of such peaks 
are sensitive to the value of the cosmological parameters, in particular to fltot, and the 
scalar spectral index n^. 

So the CMB power spectrum provides information on combinations of fundamental cos- 
mological parameters. The physical processes responsible for the anisotropics are well un- 
derstood and allows a complete prediction on the shape of the CMB power spectrum for a 
given cosmological model, see Refs. |H O [HI Ej. 
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Figure 3: The cosmic age to in units Hq^ versus flmo- The blue sohd curve describes a 
Universe with the cosmological constant, the green dashed curve corresponds to an open 
Universe without a cosmological constant, the black dotted line gives the lower bound of the 
cosmic age from the globular cluster. The red stripe gives the constrains on the age of the 
Universe from WMAP 5-year (assuming the ACDM). 
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Figure 4: The CMB power spectrum £ {1 +^)C^/27r versus the multiple £ and the angu- 
lar size 9 The curve represent the theoretical prediction of the power spectrum whereas 
the black points are the WMAP 5-year data. Taken from the webpage of WMAP: 
http://map.gsfc.nasa.gov/, 
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Generally, before recombination, the baryons in the Universe were locked to the photons 
(via Thompson scattering) of the CMB, and the photon pressure interacting against gravita- 
tional instability, due to matter, produced sound waves in the plasma. After recombination, 
the baryons and the photons separated, but the effects of the acoustic oscillations remained 
imprinted in the spatial structures of the baryons and eventually on the dark matter. The 
physical length scale of the acoustic oscillations depends on the sound horizon of the Uni- 
verse at the epoch of recombination, where the sound horizon is the comoving distance that 
a sound wave can travel before recombination and depends only on the baryon and matter 
densities. The relative heights of the acoustic peaks in the CMB anisotropy power spectrum 
measure these densities with excellent accuracy. 

The CMB anisotropics can be thought of as fluctuations in temperature 5T/T around 
the mean black body temperature T ~ 2.726K of the cosmological radiation. During the 
radiation dominated era the equation describing the effective temperature fluctuations 5T 
of the CMB has the form of a harmonic oscillator sourced by a function of the gravitational 
potential (-0), the amplitude of these oscillations are modulated by the photon-baryon sound 
speed Cg — 1/3(1 -I- ipi,/ Pr)- Therefore when photons decouple from baryons their energy 
carry an imprint of such oscillations. The characteristic frequencies of these oscillations are 
fixed by the size of the sound horizon at the decoupling Xgh = J Cgdr. Therefore we have 
a series of compressions and rarefactions at scales km^sh = mi:. Today such scales appear 
at angles that are multiple integers of the angular size of the sound speed horizon at the 
decoupling 9hs = Xsh/ Dk{tiss)i where Dk{tiss) is the distance to the last scattering surface 
for a spacetime with curvature K. As a consequence, the position of the acoustic peaks in the 
power spectrum depends on the geometry of the Universe. For a flat cosmology the peaks will 
appear at the multipoles = mish- However the acoustic oscillations are perturbed by the 
evolution of the gravitational potential which shift the position of the peaks by an amount 
that depends on the cosmological parameters that are relevant before recombination. Then, 
it is easy to understand that the position of the peaks and their amplitudes give information 
on cosmological parameters. 

In Fig. (|4]) we show the predicted CMB temperature anisotropics £{! + £) Ci/2Tr versus 
the multiple £ together with the WMAP 5-year data, see Ref. [3]. 

5.4 Baryon acoustic oscillations 

The acoustic peaks in the CMB power spectrum are predicted to be present also in the 
clusters of galaxies as a series of weak fluctuations (as a function of scale), usually called 
wiggles, see Ref. [1 H [101 El [H] . 

Measuring the acoustic peaks in the galaxy power spectrum at high precision, matching 
those already measured in the CMB power spectrum, would provide a spectacular confirma- 
tion of the standard cosmological model in which mass over densities grow from the initial 
fluctuations present in the CMB. In Fig. ([5]) are plotted the power spectra divided by a 
smooth power spectrum (without baryon oscillations). 

The position of the peaks and troughs in Fourier space can be calculated form linear 
physics and they are considered as a standard ruler, see Refs. |131I14| . Therefore an analysis 
of the galaxy power spectrum in a redshift survey containing acoustic peaks can be used to 
measure cosmological parameters. BAO have been detected at low redshift in the 2dFGRS 
and SDSS galaxy sample, see Refs. ]W[ UM [T7] . 

What is a standard ruler? 

The comoving size of an object at some redshift along the line of sight (A||) and the 
transverse (A^) direction are related to the observed sizes Az and A9 of the object by the 
Hubble parameter H and angular diameter distance Da- When the true scales, A|| and Aj^, 
are known, measurements of the observed dimensions, Az (deep in space) and A9 (wide in 
angle) , give estimates of the Hubble parameter and the angular diameter distance: these are 
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Figure 5: Ratio of the full observed power spectra to smoothed spectra to show the oscil- 
lating BAO coniponent. The power spectra are calculated from (a) the combined SDSS and 
2dFGRS main galaxies, (b) the SDSS DR5 LRG sample and (c) the combined samples (solid 
symbols with la) The solid line in each panel is the fit to the data. From Refs. [?] . 
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the standard rulers. 

In turns the Hubble parameter and the angular diameter distance (integral over the 
redshift of the inverse of Hubble parameter) are related to the matter content of the Universe, 
and information on the cosmological parameters can be extracted. 

5.5 Large scale structures 

The conventional paradigm for the formation of structures in the Universe is based on the 
growth of small perturbations due to the gravitational instabilities. In this picture, some 
mechanism is required to generate small perturbations in energy density in the very early 
phase of the Universe. The central quantity that it is needed to describe the growth of 
structures is the density contrast, defined as: 



which describes the change in energy density p{t,x) compared to the background p{t). 

Since one is often interested in statistical behaviour of structures in the Universe, it is 
conventional to assume that S and other related quantities are elements of an ensemble. The 
most accredited model of structure formation suggests that the initial density perturbations 
in the early Universe can be represented as a Gaussian random variable and a given initial 
power spectrum. The last quantity is defined through the relation: 



where S{k) is the Fourier transform of d{t,x) and (...) means averaging over the ensemble. 
Though gravitational clustering will make the density contrast non Gaussian at late times, 
the power spectrum (and the correlation function) remains the most important quantity that 
needs to be studied in structure formation. 

When the density contrast is small, its evolution can be studied by linear perturbations 
theory and each Fourier mode Sk {t) will grow independently. 

We can say that the large scale structures of the Universe start to grow after matter- 
radiation equality (teq or equivalently Oeg) when the pressure support due to photons be- 
comes smaller than the gravitational attraction due to the non relativistic matter compo- 
nent. Since non-relativistic matter has a negligible pressure relative to its energy density, 
the gravitational attraction becomes stronger and objects in the Universe start to form. The 
matter-radiation equality fixes the position of the peak in the matter power spectrum, see 
Fig. ([HI); the wavenumber keq characterizes the border of the "large scale" and "small scale" 
modes and it corresponds to the one that entered the Hubble radius at the radiation-matter 
equality, i.e. k^q = a (aeq) H (aeq). 

The observations of large scale structures such as galaxy clustering provides another 
independent observational test for the existence of dark energy. In Fig. (O we plot the 
matter power spectrum for two different flat cosmologies: ACDM where fimo = 0.3 and for 
CDM where il.mo = 1, evaluated with GAME code, see Ref. [30] 

For instance the peak position of the matter power spectrum Pm (k) is shifted toward 
larger scales (smaller k) in the presence of the dark energy. Hence the scale of the peak 
position can be used as a probe of the dark energy. 

In Fig. ([7]) are shown the galaxy power spectra for the Luminous Red Galaxy (LRG) 
and main galaxy samples of the SDSS, see Ref. [31]. As it can be seen from the figure 
the position of the peaks are around 0.01/iMpc < k < 0.02/iMpc confirming that ACDM is 
favored over the GDM model. It is worth remembering that the galaxy power spectra alone 
do not provide tight bounds on the dark energy density parameter floEo but the key point 
is that the observations of large scale structure are consistent with the existence of the dark 
energy, see Ref. |32) . 



S {t, x) 



p{t,x) ~ p{t) 
Pit) 



(55) 




(56) 
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Figure 6: The matter power spectrum in the flat ACDM model with — 0.289 and in the 
CDM model with fimo = 1- The power spectra are evaluated with CAMB code. 




Figure 7: The measured power spectra with error bars for the full luminous red galaxies 
(LRG) and the main galaxy sample of the 2dF survey. The solid curves show the theoretical 
prediction for the ACDM in the linear perturbation theory. The dashed curves include the 
non- linear correction to the matter spectrum. From Ref. |31] . 
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6 A bit of dark energy in theory 



6.1 Cosmological constant A 

The cosmological constant A is the simplest candidate for the dark energy because it is 
characterized by a constant energy density in both time and space and its parameter of 
equation of state w is constant and equal to —1. It comes also as quite natural term to be 
added in the Lagrangian, it appears in fact as a constant and linear term to R (Ricci scalar) . 

From the point of view of particle physics the cosmological constant is considered to be 
the energy density of the vacuum. 

As previously said the Einstein tensor G^i/ satisfies the Bianchi Identities V ^G^^ — 
and the energy conservation is automatically verified V^T'^" = 0. Since the metric tensor 
g^i, is constant with respect the covariant derivative there is a freedom to add terms in 
the Einstein equation proportional to g^^, without violating the total invariance. Then the 
Einstein equation given in Eq. ([I} can be written as: 

-Rm ^ ^5m"^-R + = SttGT^j. (57) 

where A is a constant. In a FLRW metric the modified Einstein equation gives: 

^ o SttG K a 

^=-^(p + 3,) + ^. (59) 
a 3 6 

It is clear to see from Eq. (|59|) that A contributes positively acting, at the background level, 
it acts as a repulsive force against gravity. 

6.2 The first problem: fine tuning 

Let us consider again the Hubble parameter: 

^ . SttG K a 

In order to have a phase of accelerated expansion today we require cosmological constant A 
to be of the order of the square of the Hubble parameter today, that is: 

A w = (2.13/1 X lO-^^GeV)^. (61) 

In terms of energy density this corresponds to: 

PA = —-^ ~ 10-4^GeV^ (62) 
an 

where Wpi = = 1.22 x lO^^GeV is the Planck mass. As previously said, the cosmological 
constant is associated with the vacuum energy of an empty space. The vacuum energy 
density is evaluated by the sum of the zero point energies of a quantum fields of mass m: 



Pvac 



(2tt) 2 



-y/k^+m^ (63) 



where k is the momentum. Of course such an integral shows a ultraviolet divergence: pyac ~ 
k^. Our belief is that field theory is valid up to a certain maximum value of frequency 
kmaxi moreover we expect General Relativity to hold just before the Planck scale: mpi = 
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1.22 X lO^^GeV. Hence if we choose k^ax to be the Planck scale then the vacuum energy 
density will be: 

Pvac WGeY\ (64) 

The ratio between the theoretical and the observed value of the cosmological constant is: 

Pvac.ob 10-47GeV'^ 

The theoretical value then is about 121 order of magnitude larger than the observed value. 

This ratio can be reduced if we assume a different cut-off scale k„iax- We expect that 
in supersymmetry theories for a certain number of bosonic degrees of freedom there will be 
exactly the same number of fermionic degrees of freedom so that the total contribution to 
the zero point energy will vanish. If the supersymmetry is broken, as it certainly is nowadays 
because we do not live in a supersymmetric Universe, then the vacuum energy will be non- 
zero. It is known that supersymmetry is broken at scale: Msusy ~ lO^GeV which will 
reduce the ratio approximately by 60 orders of magnitude but it is still not able to bring it 
down to the observed one. There have been a huge number of attempts in order to solve 
the fine tuning problem involving, QCD, supersymmetry, string theory, etc. a partial list is 
Refs. [33-[5n|. Anyway, none of them seems still convincing and moreover self satisfactory. 



6.3 The second problem: the coincidence 

The second problem of the dark energy is related again to its value today. Aside from its 
origin and the fine tuning shown in the previous subsection, its value today is comparable 
with the present matter energy density, in principle these two quantities should be unrelated. 
In other words, the value fl^^ is comparable (to a factor two or three) to flmo: for iio obvious 
reason. 

Let us consider the energy density for a general fluid, that is: 

p{t)^poa-'^'+^->^Po{l + zf'+-K (66) 

The time at which the matter energy density {wm — 0) coincides with the cosmological 
constant {wa = — 1) is given by: 

and for il^ ~ 0.7 and ~ 0.3 (best fit values) then Zc ~ 0.3 (remember that z — means 
today). In Fig. ([5]) is shown JIa as a function of the scale factor a for the ACDM model. As it 
can be seen, at early time the cosmological constant would be negligible, while at later times 
the density of matter will be zero and the Universe will be empty. Strangely (or hopefully) 
we are living in that short era, cosmologically speaking, when both dark matter and dark 
energy are comparable in magnitude. This is the coincidence problem. 

This is not a characteristic of the cosmological constant A only but of all the dark energy 
models. 

As we will show in the next section there have also been here attempts to solve, at least 
partially, the coincidence problem. 



6.4 Scalar field models of dark energy 

The cosmological constant considered so far corresponds to a perfect fluid with an equation 
of state parameter w = —I constant in time and space. This seems to be a good candidate 
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because, as we have seen so far, observations give a value of w close to —1. However, all the 
observations say little about the variation on time of such a parameter, so we cannot take 
for granted either that the cosmological constant is the true model or that the parameter 
of equation of state has always been constant. We can think to extend our horizon and to 
consider models in which the parameter of equation of state is a function of time. 

The scalar fields is not a new concept to physics, it comes quite naturally in particle 
physics and it can be also a good candidate for dark energy. So far an ample of dark 
energy models have been proposed: quintessence, K-essence, phantom, tachyonic, coupled 
dark matter and dark energy models, just to name the few. 

Here we are not interested on reporting the details of all the models present in the 
literature or all the details of a particular model. We are more interested on the impact 
on phenomenological cosmology or, so to say, we are more interested on the question: what 
should we do once we have a new model of dark energy? We briefly discuss the different 
main models while giving a list of papers where more details can be found. 



6.4.1 Quintessence 

Caldwell et ai, see Ref. [SF*, named quintessence a canonical scalar field with a potential 
V {(j>) responsible for the late time acceleration. The quintessence model is described by the 
action: 

S= I d^xy^ J-R + c^ +Sm (68) 



2fc2 

where the R is the Ricci scalar, — SttG, Sm is the matter actior0 and is the scalar 
field Lagrangian, defined as: 



The energy-momentum tensor of quintessence is: 



2 Si^C^) 



-9 '55''" 



-g^^d^<j,dpd, + V{d,) 



see Ref. [59 for more details. 

The energy density and the pressure are, in a FLRW background: 

1 -.^ 



p<t> 



.1 (0) 



(69) 

(70) 
(71) 



(72) 
(73) 



and the parameter of equation of state is: 

0^ + 2V (</>) 

The scalar field also satisfy the continuity equation: 

/90 + 3H {p4,+Pp) = 0. 



(74) 



(75) 



^We take into account all the forms of matter since it is present in the Universe as well. 
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We can also evaluate the Klein-Gordon equation just varying the action given in Eq. 
with respect to 0, that is: 

4> + + ^ = 0. (76) 



The Hubble parameter and its first derivative are: 



H 



-AttG 



1 



PM 



Pm + Pm + 



v{4>) 

J,2 



(77) 
(78) 



where we used the continuity equation for the matter component and the Klein-Gordon 
equation to evaluate Eq. (|78p . 

The acceleration parameter will be given by: 



a 
a 



AttG 



[PM + iPM ) + (<^2-y (</.)) 



(79) 



and in order to have a phase of accelerated expansion todajo we need cfP' < V (</>) . This 
means that we require a flat potential to accelerate the expansion. 

We can make a comparison with the inflation (where a scalar fleld is used in order to 
have an inflationary period) where two parameters are used, see Ref. ^60) : 



1 dy Y 



1 1 A^V 



IGTrG \ V d(j) ) 



SttGV 



A2 ■ 



(80) 



Inflation occurs when e ^ 1 and I77I ^ 1, these are called slow-roll conditions. If these 
conditions are satisfied we have: 



02«y(0), |0|«|3il(/.| 
We can write the equation of state parameter as: 

0' 



1 



IGnGeV ((/)) 



1^2 



(81) 



(82) 



where ^ = 4i/3H(j). The above equation also shows that the equation of state parameter for 
the quintessence can be bigger or at most equal to —1 (for positive values of potentials). In 
the slow roll limit, Eq. ([5T|l . we have p^j, V {(p) and ^ ^ 1 implying 1 + Wcf, — 2e/3. The 
deviation of from —1 is determined by the slow roll parameter e. 

There is now the problem in finding the proper potential which is able to lead to an 
accelerated expansion. For this purpose it is interesting to consider, as an example, a scalar 
field which gives a power law expansion: 



a (t) (X tP 



(83) 



and in order to have an accelerated expansion we need p > 1. Let us now consider Eq. (|78p 
with only a scalar field: 



H = -AnGct) 



J,2 



(84) 



where we have neglected matter for a sake of simplicity. From the last equation we can 
obtain a solution for the scalar field: 



H 
4^ 



1/2 



dt 



(85) 



'We have shown before as relativistic and non-relativistic matter are unable to accelerate the Universe. 
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and from Eq. (|78l) we have the potentiah 



1 + :^ • (B6) 



Solving Eq. and (gl]) we have: 



P .\nt, V{^)^Vocxp[-J^A (87) 



47rG ' \ \l P 

where Vb = p(3p - l)/87rG. 

The exponential potentials can give rise to the accelerated expansion, moreover they have 
cosmological scaling solution in which the energy density of the field follows the energy 
density of matter pm ■ 

Of course all the potential that are not steeper than the exponential potential can lead 
to an accelerated expansion; for instance a class of potentials is: 

V {(t>) = Af {p > 0) (88) 

V{(f>) = M^+P^-P exp(a0V"^pi) (89) 

where a is a positive number and M is a constant; for the first class of models the potentials 
do not possess any minimum and the field just roll down to its potential towards infinity; 
whereas for the second class of models the potentials have a minimum at which the field is 
trapped corresponding to u;^ = — 1, see Ref. |61j . 

In general it is complicated to construct a quintessence model because it is difficult to 
find a corresponding model in particle physics due to the energy scale which is too low. As 
we said, in order to lead to the acceleration today we require a rather flat potential to satisfy 
the condition e ^ 1 and jy^l <C 1 (slow roll parameters given in Eq. (I80p ). We can also give 
an upper bound to the mass field of the quintessence starting from the slow roll condition: 



m^i 1 



in V 



< 1. (90) 



The mass of the field is simply given by second derivative of the potential — d^F/d^^, 
hence the mass field of the quintessence today has to satisfy the condition: 




m^l <./^i/~i/o«10-=^'eV (91) 



which means that has to be extremely small in order to have an accelerated phase today. 
Nevertheless another problem arises because such a light scalar field could in principle (and 
easily) couple with ordinary matter; but such a coupling could lead to a time dependence 
on the constants of nature. However, there are tight constraints on the time variation of the 
constants of nature and, clearly, any viable model must also satisfy these constraints, see 
Ref. [62J. 

We should also have in mind that the quintessence models have been proposed to explain 
(or alleviate) the coincidence problem; this is because scalar field models have global attractor 
solutions which sub-dominantly track the dominant component of the cosmological fluid; 
however, it is still not clear which has to be the mechanism to lead quintessence to dominate 
at late times and fine-tuning the initial conditions is required. 

The key problem here is to find a physically motivated and self-consistent quintessence 
model able to account for a late time acceleration in the Universe, a partial list of quintessence 
models motivated by particle physics is given in Refs. |63|-|97j. 
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6.4.2 K-essence 



Quintessence is based on scalar fields which possess flat potentials to lead to a late time 
acceleration of the Universe. 

However, it is also possible that the acceleration of the Universe arises from a modification 
of the kinetic term of the scalar field. Originally this class of model was introduced to drive 
infiation at high energies, the so called K-inflation, see Refs. |981 199j : recent works extended 
this model and it was realized that it could have been also applied to the late time Universe, 
giving the name: K-essence, see Refs. [1001 1101| . 

Moreover, K-essence models were introduced as a good dark energy candidates because 
within a huge range of initial conditions the field automatically goes to a period in which 
w ~ — 1. In this scenario both initial conditions and coincidence problems are almost solved; 
however, other instabilities arise. 

The K-essense is characterized by a non-canonical kinematic term; the most general action 
for such model is 

S = J d^xV^P{^,X) (92) 

where the Lagrangian density P((j3,X) corresponds to the pressure p^; here is the scalar 
field and X = - i (V(/))^ 



From the action Eq. (|92p we can derive the energy momentum tensor: 

T^^) = -^^A^^fl ^ P^^d^^d.cj, + .g^.P (93) 

where the suffix " , X" represents the partial derivative with respect to X . The energy momen- 
tum tensor of the K-essence is the same as that of the perfect fiuid T^^, = {p + p) u^Ui, +p g^i, 
with velocity = df^4>/\/2X and pressure and energy density: 

= P{<f>,X) (94) 
P, = 2x'-^^-Pi^,X). (95) 
The equation of state of the K-essence is: 

P4, 2XP^x-P 

where we omitted the (f) and X dependence in P for simplicity. Moreover the K-essence field 
satisfies the continuity equation: 

p^ + iH{p^+p^)^Q. (97) 
We can evaluate the propagation speed of sound c^, which is defined as: 

2 ^ dp^/dX ^ Px , . 

' dp^/dX Px + 2XP,xx' ^ ' 

Equation ((96)) shows that the P plays a crucial role in determining the equation of state of the 
scalar field 0. Usually K-essence models have a Lagrangian of the form P (0, X) = / (0) L {X) 
but the choice is basically arbitrary. 

Moreover, the idea is that the field (p is considered as a low energy effective degree of 
freedom of some fundamental theory; such a model should respect some basic properties, 
for instance Lorentz invariance and causality. It has been shown that K-essence field during 
some epoch propagates with a speed larger than the speed of light, breaking clearly causality, 
see Ref. [T02] . 

At last, it might well be that the linear kinetic energy in X has a negative sign leading 
to quantum instabilities, such a field is called phantom or ghost scalar field. A partial list of 
K-essence models can be found in Refs. |103| - PT5] . 
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6.4.3 Phantom 



All the models considered so far have a characteristic equation of state parameter whose 
value cannot be less than —1. However, recent observations indicate that the value of wde 
has to be around —1 and moreover, they allow the possibility of a wue < — 1- If we require 
that the dark energy equation of state parameter to be less than —1 then from Eq. ([M)) we 
have: 

Px < 0. (99) 



The simplest model that satisfies the Eq. (|99|) is a scalar field with a negative kinetic energy: 

Pi<j>,X) = -X^Vi^) (100) 

where V {(f)) is the field potential. This particular scalar field is called phantom or ghost 
scalar field; its energy density and pressure are: 

P^ = -^0' + ^(0), Pc^ = -\^^-V{cl)). (101) 
The parameter of equation of state is: 

wa, — -. (102) 

^ 02 _ 21/(0) 

and for (f)^ /2 < V {(f>) we have < —1. 

However, there is also the problem of finding a self consistent and physically motivated 
scalar potential. In top of that, if the equation of state parameter is always less than —1 this 
will lead to a big-rip singularity. For instance, if the field potential has the form 

V{cb)^Voe-'^'^ (103) 

where fi is a constant then the equation of state will be always < —1. The big-rip 
singularity can be avoided if other field potentials are considered, for instance: 

V{(t>) = Foe-^'/"' (104) 
Vic/)) = Vb [cosh(/30mipZ)]"^ (105) 

(106) 

and these are called bell type potentials. In this case the equation of state starts from a 
value less than —1 but approaches to = — 1; this is because the phantom field climbs up 
the potential, due to its negative kinetic energy, but after a certain period of time the field 
settles at the maximum of the potential. 

Even if the dynamics of a phantom field seem acceptable at the classical level unfortu- 
nately it suffers from ultra violet quantum instabilities, this is because the energy density 
of a phantom field is not bound from below, which leads to an unstable vacuum state, see 
Refs. |114| and |115j . In order to avoid such a problem we need to have a phantom field that 
is weakly couple to any other normal field; but even if we assume that the phantom field is 
not coupled to any form of matter, they couple to gravitons which mediate the process: 

Vacuum 2(/« + 27 (107) 

(f> being a ghost and 7 a photon. In order not to have an overproduction of photons we 
required a Lorentz invariance breaking with a cut-off at the scale of few MeV, such a breaking 
still does not have a natural interpretation in the standard model, see Ref. |116| . 
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A Lagrangian with a negative kinetic term of the form P — ~X + V{(l)) not always 
leads to quantum instabilities. If higher order kinetic terms are added in the Lagrangian, 
i.e. P = —X + X^, instabilities can be avoided. However, it is difficult to use such a 
Lagrangian as a candidate for dark energy itself, this comes directly from the small value 
of the energy density of the scalar field which gives X ^ X^, though the stability cannot 
always be guaranteed. See Refs. |117| - [T^ for a list of papers on the phantom model. 



6.5 Dark energy as a fluid 
6.5.1 Perfect fluid 

So far we have considered only scalar fields as a candidate for dark energy; an entire class of 
models exists involving perfect fluids. Here we consider a single fluid model with a general 
energy momentum tensoiQ: 

= ip + p) Uf^u^ + pg^^ + [qf_,u^ + q^Uf^ + tt^i,] (108) 

where, beside the familiar symbols p, p and for the energy density, pressure and four- 
velocity, we introduced the heat flux vector and the viscous shear tensor tt^^. The terms 
appearing in the square brackets in the Eq. (llOSp are important only when the internal energy 
of the fluid is comparable to the total energy density. For a perfect fluid they are all zero. 
The covariant derivative = gives the conservation equation: 

p + 3H{p + p) = (109) 

valid at background level. From the above equation we can derive the expression for the 
energy density as a function of time under the assumption that the pressure is a function 
of the energy density p — wp. For the sake of generality, we consider w to be a function of 
time, the energy density scales as: 



po exp 



-3 / i±^da' 



ao = l 



(110) 



where po is the energy density of the fluid today. 
Eq. piOp can be written in a more familiar way: 

where 

^(a) = ^ r^da'. (112) 



At background level a perfect fluid is characterized by its equation of state parameter w. 
There is now the problem to parameterise w (a). As realized from the previous discussions, 
there are different possibilities: according to the preferred model the equation of state pa- 
rameter can be a constant (as it was for the cosmological constant) or it can be a function 
of time (as it was for the scalar field) . There is not a unique general expression for the dark 
energy parameter of equation of state if we consider the dark energy as being a fluid. The 
easiest thing to do is to Taylor expand w (a) around = 1; we have, see Refs. [2261 1229j 
(usually called CPL paramatererization) : 

2; 

w (a) ^ Wo +Wail - a) = Wq +Wa— — (113) 

1 + z 



^We apologize for being sometimes too repetitive. 
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where wq is the value of the parameter of equation of state today and Wa gives the variation 
on time of w (a). The values of the parameters are usually assumed to be — 1 because 

observations tell us that the value of the parameter of equation of state today has to be 
close to —1 and Wq ~ or small enough because we do not have any strong evidence that 
the dark energy varies so much over time. Moreover, the expression given in Eq. (jll3p is 
very practical because it also avoids the complication of having unrealistic behavior, e.g., 
w < -1. 



6.5.2 Chaplygin gas 

Another model involving a fluid is known as Chaplygin gas. Its pressure takes a simple form: 

(114) 

P 

where ^ is a positive constant. The continuity equation Eq. ([5]) gives: 



where B is also a constant. The expression of the energy density for the Chaplygin gas has 
two interesting asymptotic behaviors: 

/^\ 1/6 

p ~ when a < ( — I (116) 



A 



P 



5X1/6 



\fA when » > ( ^ ) ■ (H''') 



The Chaplygin gas behaves as pressureless dust at early times (when a is small) whereas at 
late time (when a ^ 1) it behaves as dark energy with w ^ —\ leading to an accelerated 
expansion; hence the Chaplygin gas can be also a good candidate for unified dark energy and 
dark matter models. However, observations exclude the Chaplygin gas being a good candi- 
date for dark energy. This is because the sound speed (the velocity with which perturbation 
propagates in the fluid) of such a model takes: 

dp B + Aa^ \^^^\ ^ ' 

hence it behaves as pressureless dust in the Matter Dominated Era (a ^ 1) and it approaches 
the speed of light at late times (a ^ 10. A sound speed of the order of the speed of light 
has implications on the growth of perturbations. A non zero sound speed gives a non zero 
Jeans scale (characteristic scale for the growth of perturbations); the higher the sound speed 
the larger is the Jeans scale. All the perturbation smaller of such characteristic scale cannot 
grow, in particular they should manifest an oscillation behavior. The fact that the latter are 
not observed in the matter power spectrum has led people to rule out the Chaplygin gas as 
a competitive candidate. It was also found that combining the results from CMB and SNIa, 
the Chaplygin gas is ruled out at more than the 99% confidence level, see Ref. |130J 

This situation can be partially alleviated if we consider a generalized Chaplygin gas; in 
this case the pressure becomes: 

P=-4 (119) 

with < a < 1. However, a has been constrained to be \a\ < 10^^ at 95% confidence level, 
see Ref. [131j. Then it seems that ACDM model is the prefered model. A partial list of 
generalized Chaplygin gas models is Refs. p32| - P^5] . 

^"The importance of the sound speed will be clear when we discuss the perturbation of a fluid. 
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6.5.3 Coupled dark energy 

As we said there is still the problem of explaining why the value of the dark energy density 
today is comparable with the present matter energy density. This fact may suggests that 
the dark energy is coupled to dark matter. The problem is finding a physically motivated 
coupling term. Let us recall the Einstein equation: 

G^, = SttGT^,, (120) 

we know that from the Bianchi Identity, the Einstein tensor G^^ is conserved which guaran- 
tees the conservation of the energy momentum tensor: 

V^T^, = 0. (121) 

However, the energy momentum tensor that appears on the right hand side of Eq. (|120p has 
to be considered as the sum of the energy momentum tensors of all the species present in the 
Universe, T^^, = X^i^M^'- Usually it is taken for granted that all the species are decoupled 
(but always coupled through gravity) so that the Eq. (|12ip holds independently for all the 
species: V^r^l} = 0. However interactions between different fluids, say matter and dark 
energy, might be likely to happen. Hence Eq. (|12ip (which is always valid) is the starting 
point for finding the coupled model. 

Several coupling forms between dark energy and dark matter have been discussed in 
literature. Amendola proposed, in Ref. ||150j , a quintessence scenario coupled with dark 
matter as a direct consequence of non minimally coupled theories; the coupling term has the 
form Qpm4>- 

Another interaction, which is based on the fluid description, is Hpm6, where (5 is a 
dimensionless coupling, see Ref. [HI], [Hg, [1131, [IH, [HSI, [ISSI- 

As an example, let us consider the dark matter and the dark energy, coupled together, 
with an interaction of the form: 

p„ + 3Hp,n = Hp,nS (122) 

Pde + iH [pde + Pde) - H p,n5 ■ (123) 

The coupling 5 can either be constant or time varying 
• Constant 5 

The Eqs. (|122p and (|123p can be solved easily, giving: 

Pm(a) - Pm„a"'+^ (124) 

Pde[a) = Pde„a ^ ^ + J— 

+ 3Wde 



(125) 



where we assumed Wde to be constant. The Hubble parameter then reads: 



Hi 5 + iWde 



3u;dea"^+'' + <5a"^^'+"'^ + r!<ie„a"^^'+"'''=^ (126) 



Observational constraints on the coupling constant 5 are plotted in Fig. combining 
SNIa, CMB and BAO. 

• Time varying 6 

In order to have an expression for the Hubble parameter we need to find a parame- 
terisation of 5{a). Generally speaking, we are faced with two possibilities: either we 
write down a parameterisation for 5(a) and then find from this the relation between 
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a 



log(a) 



Figure 8: Plot of the fl\ versus log a assuming flat space with rt\„ = 0.7. There are 
also indicated the Planck era, the electroweak phase transition (EW) and the Big Bang 
Nucleosynthesis (BBN). From Ref. [FT] 




Figure 9: Probability contours at la, 2a and 3a respectively combining SNLS, BAO and 
CMB data. The left panel shows observational contours in the [wx , S) plane marginalized 
over flxa, whereas the right panel shows contours in the (flxoi^) plane marginalized over 
Wx- The best-fit model parameters correspond to S = —0.03, wx = —1.02 and flxo = 0.73. 
Form Ref. [157j 
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Pmifl) and pdeia), or first give the latter and then derive a function (5(a). Here we show 
this second approach. The basic relation we start from is that the fluid densities scale 
according to the following relation, see Ref. [153] 



Pdeja) 
Pm (a) 



Pdeo f 
Pmo 



(127) 



where ^ is a constant parameter. Moreover, we approximate Wde as a constant: it is 
clear however that a more complete analysis should allow for a time-dependent equation 
of state. 



The relation given in Eq. (|127p has two useful properties: a) it includes all scaling 
solutions (defined as those with p ~ a™) and b) the functions Pm(a) and pdeio) can be 
calculated analytically. Since for uncoupled dark energy models with constant equation 



of state one has p„ 



and pde 



reduces to this case for ^ = 
deviates from the law. 



- a-3(i+«;de) it appears that the relation Eq. (fT27| 
Conversely, if ^ 7^ —iwde then the matter density 



To find a parameterisation of 5{a) we can derive with respect to time Eq. (|127p and 
substituting Eq. and Eq. we find then: 



5 (a) = -^de (a) (3wde + 



(128) 



where ^de{a) = Pde{a) / {pde{a) + Pm{a)) = ^deg/ {^deg + ^mga ^), hence we can write 
Eq. (fT28l) as: 

^(«) = ,0 (129) 



deo 



where 6q — ^fldeo (Sw^e + 1^) is the present value of the coupling. 

Using Eq. (|122p and (|123p with the Eq. p27p we can find the conservation equation 
for the effective fluid ptot = Pm + Pde, that is: 



ptot = -iH 



1 + u; 1 



Pdeo 



Ptot 



(130) 



integrating the latter we find the Hubble parameter: 



(131) 



There are three parameters {wq, ^, ^deg} to constrain, we can use 5q instead of ^ since 
they are related via (5o = ~^deo {Swde + £,)■ Data analysis using SNIa, CMB and BAO 
gives the bounds -0.4 <S < 0.1, -1.18 < Wde < -0.91 and 0.69 < Qdeo < 0.77 at the 
95% confidence level, see Ref. |157| . 



7 Dark energy: modification of gravity 

Because of the unfeasibility of solving the theoretical problems related to dark energy, cos- 
mologists have started to look at other models. An alternative approach postulates that 
General Relativity is only accurate on small scales and has to be modified on cosmological 
scales. We have seen that the contribution of the matter content is given by the energy 
momentum tensor on the right hand side of the Einstein equation, whereas the left hand 
side represents the geometry. In order to have a late time accelerated phase we included an 
extra component, called dark energy, in the Universe which is meant to add an extra term 
on the right hand side of Einstein equation. However, an effective accelerated phase can be 
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also reproduced by modifying the geometry of the Universe or by simply change the laws of 
gravity. To understand that, let us assume that the Newtonian gravity is described by: 



F ^ G^e-^/^o (132) 

where m and M represent the masses of two objects, r is their distance and vq is some 
typical fixed distance. For small value of r Eq. (|132p becomes similar to the traditional 
Newtonian force, i.e. F ~ G ™ and we are tempting to say that the last one represents 
the correct theory but if we look a little bit further, say at distances comparable to rg then 
the gravitational attraction becomes weaker as a consequence we are led to say that there 
has been an extra force which has counterbalanced the gravitational attraction. This is what 
happens in the modified gravity models. 

This review is not meant to explain in detail the different models existing nowadays but 
for the sake of completeness we will briefiy list the most studied modified gravity models, 
pointing out the main characteristics useful for the phenomenological cosmology. 



7.1 f{R) gravity 
The Einstein-Hilbert action 



S 



r 1 
2k 



R + C 



-gd X 



(133) 



is the base of General Relativity, by simply varying the latest with respect to the metric we 
can obtain the Einstein equations seen before. 

The / (i?) gravity intends to modify the action itself; in principle there is no a priori 
reason to consider a gravitational Lagrangian as a linear function of the Ricci scalar R. Let 
us consider an action of the form, see Refs. |159| 1160] : 



S= f{R) V^d^x 



(134) 



where R, in the Einstein-Hilbert action, has been replaced by a generic function / (i?) of R. 
The field equation can be obtained by varying the action (I134p with respect to the metric 

(135) 



f{R)R^, - -fiR)g^u - [V^ V, - .g^.nj/'(i?) - SnGT^, 

where □ = g^iud^d'^ is the d'Alembert operator, the prime refers to the derivative with 
respect the Ricci scalar i?, i.e. f'{R) = df{R)/dR, and T^^ is the matter energy momentum 
tensor. Taking the trace of the Eq. (|135|) we have: 



3n/'(i?) + Rf{R) - 2f{R) = 8nGT 

being T = —p + 3p the trace of the energy momentum tensor. 
Working in the flat FLRW metric, Eqs. (fTSSj) and (fTSe]) read: 



(136) 



2H + 3H^ 



1 [Rf'iR) f{R)] - SHRf'iR)] (137) 



SttG 



i(/(i?)-i?/(i?))}. 



(138) 
(139) 
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As we said f{R) theories are interesting because they lead to accelerated expansion without 
the need for dark energy. To see this we can define an efliective energy density and pressure 
as: 



1 

Peff 



Rf'{R)-f{R) --^HRf"{R) (140) 



2/'(i?) L ^ ^ ' ^ ^ 'J /'(i?) 
Peff = -j^^[ief''{R) + 2HRf'{R) + Rf'{R) + \[f{R)-Rr{R)]] (141) 

where pef / has to be positive to guarantee a positive Hubble parameter (one way to see this 
is to set in Eq. (jl37p p = 0, hence, in vacuum the curvature correction can be viewed as an 
effective fluid). 

The effective equation of state parameter is: 

R^r{R) + 2HRf"{R) + Rf"{R) + \ [f{R) - Rf{R)] 
"^eff = — — — — — ■ (142) 



Rf'iR)-fiR) -iHRf"{R) 



The effective parameter of equation of state depends on the function f{R); let us give two 
examples: first, we consider a f{R) oc i?" and the scale factor is assumed to be a generic 
power law a{t) — {t/tQ)", we obtain, see Ref. |159) : 

6n2-7n2-l 

"^^^ = - 6n^-9n + 3 ^'''^ 

for n 7^ 1 and a is given by: 

-2n2 + 3n-l , 

a = . 144 

n — 2 

According to the value of n a desired value for Wef / can be obtained. For instance if n = 2 
then Weff = — 1 and a = oo as desired. However, this model was first considered as candidate 
of infiation where / (i?) = R + aR^ which leads to an accelerated epoch but it ends as soon 
as aR^ < R, see Ref. 



The second example uses the function f{R) = R - where fi IS EL generic 

parameter, see Ref. |159| : in this case, assuming the scale factor to be a power law, we have: 



2 (n + 2) 
3(2n+l)(n+ 1)' 



^e// = -l+ ,,.„ : (145) 



For instance, if n = 1 then we have Weff = —2/3, see Ref. [160J. However, also this class 
of models are not good candidates for the late time acceleration for several reason: first of 
all they are excluded by solar system constraints (this is because /" < 0), see Ref. [1611 
11621 11631 11641 1165j : then, in the Ref. 1167j was shown that this class of models lead to 
a different growth of the scale factor during matter dominated era, hence inconsistent with 
observations. 

It is also worth mentioning that some of the problems listed above can be solved if the 
Palatini formalism of / (i?) gravity is used; in such a formalism the independent variable 
is not anymore the metric tensor g^j, but the connection FJ^^. As a consequence the Ricci 
tensor calculated in the Palatini formalism i?^^ (F) will differ from the one calculated in 
terms of the metric tensor R^i, (g). Even though all the stabilities problems are solved at the 
background level, there is still the problem of the evolution of matter density perturbations 
which may manifest a violent growth or a series of damped oscillations being incompatible 
with observations, see Ref. |168) . 

Good review of f{R) theory can be found in Ref. [169J and |170| . 
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7.2 Brane world: DGP 



It might also well be that we live in an Universe with more than four dimensions (three 
spatial dimensions and one temporal dimension); although the other extra dimensions are 
compactified on some manifold in order to obtain 4-dimensional gravity theories. One of 
the best studied examples is the Dvali-Gabadadze-Parroti brane- world model (DGP), see 
Ref. |171j . In this model the 4-dimensional brane is embedded in a 5-dimensional Minkowski 
bulk with an infinitely large extra dimension. Matter is confined in the 4-dimensional brane 
and only gravity can propagate in the 5-dimensional bulk. Gravity leaks off the 4-dimensional 
brane into the 5-dimensional bulk only on large scales, say A > r,,. On small scales, A < 
gravity is effectively bound to the brane and 4-dimensional dynamic is recovered with very 
good approximation. Then the gravity leakage leads to the observed late-time acceleration 
of the expansion of the Universe. In other words the acceleration is not due to the presence 
of a dark energy component but rather to the weakening of gravity. 
The action, for this model, is: 




S - -— ^ / d^yy^i?5 - / d'*a;V-/ii?4 + / dV-Z^An + ^ghy (146) 



where gap is the metric in the bulk and ft.^^ is the induced metric on the brane; Cm is 
the matter Lagrangian confined, as we said, on the brane. The term 5ghy is the Gibbons- 
Hawking- YorlOboundary term which needs to be added all the time the space time manifold 
has a boundarjo. 

The characteristic scale we said before is given by the ratio of the Planck mass A/pi and 
its counter part in the 5-dimensional bulk, that is: 

The action (|146p gives the following modified Hubble equation: 

H^-^H^ (148) 

where e ~ ±1, this because the DGP model gives two different solutions, we will consider 
here only the self-accelerating one: e = 1; p is the total energy density of the fluid confined 
on the brane and it satisfies the standard conservation equatiorP^: 

p + iH{p + p)=Q. (149) 

Eq. (|148p can be written, assuming only matter present on the brane and zero curvature, as: 



where r2,„ is the matter density parameter and 



(150) 



1 

'0 



Setting z = in Eq. p3D)) we get: 



^r. ^ y . (152) 



^^The necessity of such a term was first pointed out by York and later it was refined by Gibbons and 
Hawking. 

^■^ Adding this extra term the variational principle is well-defined. 

^•^Matter is confined on the brane and hence there is no flow of matter outside the brane. 
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Comparing Eq. (jl48p to the normal Friedmann equation with an additional dark energy 
component, we see that we can move the crossover term to the right hand side and think 
of it as a dark energy contribution with pde — iH / {SirGrc). Looking at the conservation 
equation (|149p we find that it is solved if the effective dark energy has an equation of state 
given by 

l + ^DE--3§2. (153) 

However, things are not as nice as they look. In DGP model there has to be a sort of fine- 
tuning as welO: we need to quantify the cross-over scale Tc- If we consider the equation 
Eq. p48p we realize that in a CDM domination situation characterized by p cx the 
Universe will approach the de Sitter solution; 

H^H^^-. (154) 

This means we can have an accelerated phase at late time without invoking dark energy 
but the price to pay, because we require the accelerated phase to happen now, is that the 
cross-over scale has to correspond to the present Hubble radius today, i.e. Tc ^ Hq. With 
this fine tuning the DGP model seems to fit the SNe observations; however, when we add 
other observations such as BAO and CMB shift the DGP model is ruled out, see Ref. |172] . 

However, the DGP model itself suffers from several stability problems (namely ghost), see 
Ref. |173) - [T75] : moreover it is also about to be excluded from observations, especially from 
solar system tests, see Ref. [176j and |,177j . These difficulties have led people to extend the 
DGP model to higher dimensions, called cascading DGP, see Ref. |178| and |179| . Recently, 
people have developed a covariant cascading gravity in 61? brane-world model to study the 
cosmological solutions, see Ref. |182| . 

Another way to keep the DGP model ghost-free is to assume a light scalar field which 
is coupled to metric tensor; furthermore, it is fundamental to have only second order time 
derivatives in the field equations otherwise this would rise to extra degree of freedom which 
are likely to give rise to a ghost state. However, these theories might have superluminal 
excitations, see Ref. [183] and [184'|. 

A review of brane-world cosmology can be found in Ref. [185) . 



7.3 Neither of the two 

Instead of modifying the right hand side of Einstein equation (adding dark energy) or the 
left hand side (modifying gravity) to explain the observed late time acceleration, a third 
possibility is to drop the assumption that the Universe is homogeneous an large scales. 

The assumption that the Universe is homogeneous and isotropic is known to be violated 
at late times when structures formation starts to be non-linear. Therefore, it is important 
first to quantify the effect of the breakdown of linearity before invoking new physics. 

The main idea comes from the fact that the Friedmann equation is derived from the 
Einstein equation G^i, — SttGT^^ which is non-linear in the metric g^,y (the Einstein tensor 
G^jy contains quadratic and quartic terms in the metric). The Universe expansion and the 
scale factor are supposed to describe the average dynamics of the Universe, hence they should 
be evaluated from the Einstein equation averaged over large scales: 

(G^,) =87rG(T^,). (155) 

However, the approach of the FLRW model is to define first the average metric (5^,^) — 
diag(— a^, a^, a^, a^), an average density and pressure and then to compute the Einstein 
equation, that is: 

G^. i{g^.)) = 87rGr^^ ((p), (p)) . (156) 

^*Problem that we tried to solve in the ACDM scenario. 
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This equation leads to the Friedmann equation, but clearly Eq. (|155p and Eq. (|156p do not 
give the same results because of the non linear structure of G^i^. Moreover G^^, contains 
time derivatives and spatial derivatives of g^^; hence, it is possible that g^v is very close to 
{9^.v) while at the same time Gfn, and {G^u) may be significantly different. Fortunately, this 
effect can be neglected before the time of equality of matter and radiation, because it is clear 
that the Universe is nearly homogeneous before that time, but departure from homogeneity 
should be taken into account when structures start to form. The majority of cosmologists 
believe that the standard approach provides a good approximation, relying especially on the 
N-body simulations. A few experts challenge this picture and try to compute the proper 
averaging, see Refs. [1861 1187| . 

This idea is of course the most appealing and economical one that one might think of, 
but there are doubts that it could work. The problem is that technically it is impossible to 
solve analytically the Einstein equation (if special symmetries are considered) and moreover 
N-body simulation is still based on Newtonian gravity. 

However, the idea of considering the effect of inhomogeneities has led people to consider a 
different metric instead of the usual FLRW one. One example is the Lemaitre-Tolman-Bondi 
(LTB) model. In this model the Milky Way has to be near the center of a large underdense 
region, called void, in order to be consistent with SN la observations and CMB, see Ref. |180j 
and |181j- Of course in this case we need to violate the Copernican principle requiring that 
we live a privileged location, see Refs |188[ 1189) . 



8 Dark energy or modify gravity? 

There is now strong observational evidence that the expansion of the Universe is accelerating. 
The standard explanation invokes an unknown "dark energy" component. But such scenarios 
are faced with serious theoretical problems, which has led to increased interest in models 
where instead General Relativity is modified in a way that leads to the observed accelerated 
expansion. The question then arises whether the two scenarios can be distinguished. Here we 
show that this may not be so easy. We demonstrate explicitly that a generalized dark energy 
model can match the Dvali-Gabadadze-Porrati model and reproduce the 3 + 1 dimensional 
metric perturbations. Cosmological observations are then unable to distinguish the two cases. 



8.1 Background quantities: is this all we need? 

It is straightforward to explain the observed accelerated expansion of the late-time Universe, 
within the framework of FLRW cosmology by simply introducing a cosmological constant or 
a more general (dynamic) dark energy component. 

An alternative approach postulates that General Relativity is only accurate up to certain 
scales and has to be modified on cosmological distances. This in turn leads to the observed 
late-time acceleration of the expansion of the Universe, see Ref. [1911 1192|, I193| . 

One important question is whether a scenario where the gravity is modified can be dis- 
tinguished from one invoking an invisible dark energy component. It is well known that any 
expansion history (as parameterised by the Hubble parameter H{t)) can be generated by 
choosing a suitable equation of state for the dark energy (parameterised by the equation of 
state parameter w — p/ p oi the dark energy): 



2 {1 + z)H{z)H'{z)~ H{zy 

H{zY - Hln^ (1 + zf 



where the prime refers to the derivative with respect to the redshift. Let us illustrate this 
explicitly for the DGP model, for which the Hubble parameter evolves as Eq. (|148p ; then the 
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effective dark energy equation of state is given by: 

1 H'{z'\ 

w{z) = -l + -{l + z)^. (158) 

Consequently, it is impossible to rule out dark energy based on measurements of the cosmic 
expansion history (e.g. SN- la data). 

Lately there have been claims that if we are able to measure the growth rate of structures 
then we can rule out different models. This because different theories predict different growth 
of structure, then fixing the parameter of equation of state w (effective or not) we can predict 
the evolution of the growth rate; if we observe a different growth then the theory is wrong 
and can be ruled out. 

8.2 The role of perturbations 

Here we want to evaluate the growth of matter perturbations in two different scenarios: 

• General relativity and dark energy 

• Modified gravity (without dark energy) 

in order to verify the previous statement that at perturbations level different models, with 
the same equation of state parameter, can be distinguished. 

In order to compare models of dark energy with observations like the Cosmic Microwave 
Background (CMB) and large scale structure (LSS), it is important to study the evolution 
of density perturbations in a Universe containing dark energy. 

As mentioned, the conventional paradigm for the formation of structures in the Universe 
is based on the growth of small perturbations due to the gravitational instabilities. The 
central quantity is the density contrast 5 {t, x), defined in Eq. ([55)) . 

The growth of perturbations depends on the characteristics of the fluid itself, namely 
pressure, pressure perturbation and perhaps anisotropic stress. For ordinary matter this 
quantities are known, so the growth of the density contrast can be easily evaluated. However, 
for dark energy these parameters are a priori unknown and need to be measured. Usually it 
is taken for granted that dark energy does not cluster on small scales because it is assumed 
that its sound speed is too high to allow perturbations to gravitationally collapse. For some 
specific models, i.e. scalar field models, this is a good approximation but for a general dark 
energy fluid this may not be the case and perturbations in energy density have to be studied. 

The dynamic evolution of dark energy can affect a number of observables, including 
the spectrum and the growth of large scale structures, weak gravitational lensing, SNIa 
apparent luminosity and CMB anisotropics. One particular manifestation occurs in the late 
time Integrated Sachs- Wolfe (ISW) effect, which measures the evolution of the gravitational 
potential as the Universe enters a phase of dark energy domination. This effect is only 
significant on large scales (low multiples), since small scale fluctuations in the gravitational 
potential smooth out along the line of sight. And it is only significant at late time since 
potentials evolve the same as the background during matter domination. The ISW effect 
has been detected in cross-correlations between CMB temperature anisotropics and surveys 
of large scale structure, see Ref. |195| . 

It is imperative to deeply study the perturbations equations in order to compare different 
models with observations, as they clearly depend on them. 

Here we provide the perturbation equations in a dark energy dominated Universe starting 
with a barotropic fluid and we then extend to a more general fluid. We will concentrate then 
only to dark energy perturbations. However, also at perturbations level, dark energy suffers 
of serious stability problems. We face these difficulties and we show the attempts made in 
order to encompass these stabilities. 



35 



For simplicity, we consider a fiat Universe containing only (cold dark) matter and dark 
energy, so that the Hubble parameter is given by 



/ \ nn \ 



where w{a) is given in Eq. (I112[) . 

We will consider linear perturbations about a spatially-flat background model, defined 
by the line of element: 

ds^ = [- (1 + 2A) drf + 2Bidr]dx' + ((1 + 2Hl) % + '^Hnj) dx.dx^] (160) 

where A is the scalar potential; Bi a vector shift; Hl is the scalar perturbation to the spatial 
curvature; HJ^ is the trace-free distortion to the spatial metric, see Ref. [196| for more details. 

We will assume that the Universe is filled with perfect fluids only, so that the energy 
momentum tensor takes the simple form: 

T'''' = {p + p)u''u'' +p gf" (161) 

where p and p are the density and the pressure of the fluid respectively and is the four- 
velocity. 

The components of the perturbed energy momentum tensor can be written as: 

To" = ~{p + Sp) (162) 
T° = {p + p)iv,-B,) (163) 

n = {p + p)v^ (164) 
= {p + Sp)S'^+pIl'j. (165) 

Here p and p are the energy density and pressure of the homogeneous and isotropic back- 
ground Universe, Sp is the density perturbation, dp is the pressure perturbation, is the 
vector velocity and H* is the anisotropic stress perturbation tensor which represents the 
traceless component of the Tj. So far the treatment of the matter and metric is fully general 
and applies to any form of matter and metric. 

As it was at background level we have to choose now a preferred metric. Here, we 
want to investigate only the scalar modes of the perturbation equations. We choose the 
Newtonian gauge (also known as the longitudinal gauge) which is very simple for scalar 
perturbations because they are characterized by two scalar potentials ijj and 0; the metric 
Eq. p60p becomes: 

ds^ = [- (1 + 2^j) dr]^ + (1 - 2(/)) dx^dx'] (166) 

where we have set the shift vector Bi — and H!^ = 0. The advantage of using the Newtonian 
gauge is that the metric tensor g^^, is diagonal and this simplifies the calculations. Moreover, 
as we will see later, it is very useful to use the Newtonian gauge for the observational tests. 

In the conformal Newtonian gauge, the first-order perturbed Einstein equations give, see 
Ref. |T97j for more details: 
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where a is related to Ilj through: 

(p + p) a = - (k,k, - 1,5,, j n;-. (171) 

The energy-momentum tensor components in the Newtonian gauge become: 

Tl = -ip + Sp) (172) 
ikiT^ = -ihT^^{p + p)9 (173) 

t; = (p + ,5p)5;.+pn;. (i74) 

where we have defined the variable 9 — ikjV^ which represents the divergence of the velocity 
field. 

The perturbations equations are obtained taking the covariant derivative of the perturbed 
energy momentum tensor, i.e. T^^^ = 0. We have: 

6 = - {1 + w) (e - 34)^ - 3^ " w6^ for i/ = (175) 

= --(l-3w)9 —9 + k^-^^^ + k^iP-k^cr for i/ = i. (176) 

' 1+w l+w ^ ' 

The equations above are valid for any fiuid; for non-relativistic matter (for instance for cold 
dark matter) we need simply to set w — 0, dp = and a = 0. 



S = -9 + 3^ (177) 

9 = --9 + e^p. (178) 
a 

The above equations depend on the evolution of the gravitational potentials which depend 
on the evolution of the perturbations of all the species in the Universe (if the cluster at all!) 
just because the right hand side of Eqs. (|167p - (|170p have to be considered as the sum over 
all the species. As we mentioned at background level, any fiuid is characterized by its own 
equation of state parameter w but at first order perturbation level we need to specify also 
6p and a in order to close the system. 

Here we want to focus on dark energy perturbations because, as we will see, they play 
an important role on the evolution of matter perturbations, hence on the evolution of large 
scale structures. 

The equation of state parameter can be inferred by measuring the Hubble expansion 
whereas dp and a depend on the intrinsic characteristic of the fluid and the choice of these 
two functions is almost arbitrary because we have not been able to detect dark energy in the 
laboratories yet. 

However, particular attention has to be paid to the pressure perturbation Sp because it 
is responsible for the growth of perturbations. 

The pressure perturbation of a particular component is related to its velocity dispersion 
w; if the pressure distribution of this component is higher than the gravitational attraction 
then the density perturbations can not growth. Only if we lower the pressure perturbation 
then the amplitude of the density contrast can increase. This is what happens for (dark) 
matter: where the pressure perturbation is basically zero (it is often said that matter is 
pressureless) ; in the last case the gravitational attraction is always bigger than the pressure 
support and the perturbations can grow at any time and any scale. 

The velocity dispersion (or equivalently pressure perturbation) will define a characteristic 
scale called the Jeans scale: A,/ — 1/kj. In brief, all the perturbations smaller than this 
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typical scale cannot grow just because the pressure support is bigger than the gravitational 
attraction; if a particular perturbation is bigger than the Jeans scale then it is free to grow 
because it does not feel the pressure support inside this region. Hence we need to characterize 
the pressure perturbation for any clustering fluid present in the Universe. 

However, the discussion is here a bit more complicated since the Jeans scale is in comoving 
coordinates and hence it grows on time; it may happen that at some time ti a perturbation 
with a scale A (or equivalently k ~ 1/A, anyway it is a fixed number) was bigger than the 
Jeans scale but at some t2 > ti the Jeans scale has become bigger than A and this particular 
perturbation it is said that it has entered inside the sound horizon and the perturbation at 
this corresponding scale has stopped growing. 

Usually all the fluids in the Universe are considered to be adiabatic (the definition of 
adiabatic fluid will be clear in the next section), so the pressure perturbation takes a simple 
form: 



where we have introduced a new quantity — p/p, called adiabatic sound speed. For 
regular matter and radiation the adiabatic sound speed is constant: — and — f/3, 
respectively. We also pretend dark energy to be adiabatic and hence to have the same 
expression for the pressure perturbation but soon we realize that if we assume the equation 
of state parameter for dark energy to be constant and negative (negative in order to lead to 
an acceleration phase) we automatically have a negative sound speed with a consequently 
exponential growth of the perturbations. Moreover, if w is a time dependent function close 
to —1, the second term in Eq. p79p will diverge. 

The latter ones force us to go beyond the assumption of adiabatic fluid. In the next section 
we will also consider a different attempt in order to stabilize the pressure perturbation. 

We should also consider the anisotropic stress; it will be a bit bias to assume it to be zero. 
The choice of a is even more difficult as it might rise from a pure geometrical configuration or 
from a real internal degree of freedom. A familiar class of models represents the anisotropic 
stress as a viscosity term of the density perturbations, hence it should be clearly connected 
to the the velocity perturbations and its first derivative 

Later we will see the importance of the anisotropic stress and how it may influence the 
growth of matter inhomogeneities. 

8.2.1 How to deal with crossing models 

We saw before that in the Newtonian gauge we have three parameters to measure in order 
to close the system and to solve the equations of perturbations: w, Sp and a. 

Current limits on the equation of state parameter w — p/p of dark energy seem to 
indicate that p —p. see Refs. |198| I199| . sometimes even that p < ~p. Ref. |200| . often 
called phantom energy, see Ref. |201| . Although there is no problem considering w < —1 
for the background evolution, there are apparent divergencies appearing in the perturbations 
when a model tries to cross the "phantom divide" w = — 1, see Ref. [202 . Even though this 
region may be unphysical at the quantum level, Ref. |2031I204] . it is still important to be able 
to probe it, not least to test for alternative theories of gravity or higher dimensional models 
which can give rise to an effective phantom energy, see Refs. j2051 12061 12071 [^TT[ I212| . 

However, at the level of cosmological first order perturbation theory, there is no funda- 
mental limitation that prevents an effective fluid from crossing the phantom divide. 

As w ^ — I the terms in Eqs. (|I75p and (II76P containing 1/(1 + w) will generally diverge. 
This can be avoided by replacing with a new variable V defined via V = [1 + w)6. This 
corresponds to rewriting the Qi component of the energy momentum tensor as ikjTg = pV 
which avoids problems if Tq ^0 when p — —p. Replacing the time derivatives by a derivative 




(179) 
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with respect to the scale factor a (denoted by a prime), we obtain: 



S' = 3{l + w)i;' -^-3-(^-ws] (180) 
iia^ a \ p J 

v ^ -a-3-)- + ^f + (1+-)^^ (181) 

where we are assuming (at the moment) the anisotropic stress a to be zero, implying (f) = ip, 
see Eq. (|T70)) . 

In this form everything looks perfectly finite even at w — —1. But we still need to give 
an expression for the pressure perturbations. 

In addition to the fluid perturbation equations we need to add the equation for the 
gravitational potential ip. If there are several fluids present, then the evolution of each of 
them will be governed by their own set of equations for their matter variables {6i,Vi}, linked 
by a common ip (which receives contributions from all the fluids j^. 

The evolution of a fluid is therefore governed by Eqs. (|180l) and pSip . supplemented by 
a prescription for the internal physics (given by Sp) and the external physics through ip and 
ip', by Eq. (fW)l . 

We start with a barotropic fluid highlighting the problems; we extend the discussion to 
a more general case trying to overcome these difficulties. 

• Barotropic fluid 

We define a fluid to be barotropic if the pressure p depends only strictly on the energy 
density p: p = p{p). These fluids have only adiabatic perturbations, so that they are 
often called adiabatic. We can write their pressure as 



dp 

pip) = pip + Sp) = p{p) + -j- 

ap 



Sp + 0{{5p)'). (182) 



Here p{p) — p is the pressure of the isotropic and homogeneous part of the fluid. 
Introducing iV = In a as a new time variable, we can rewrite the background energy 
conservation equation, p = — 3'H(p + p) in terms of w, 

dw dp dw dw 

-TTT = -TTT-;— = —J 1 + ■ (loo) 

dN dN dp ^ dp ^ ' 

We see that the rate of change of w slows down as w — > —1, and w — ~1 \& not 
reached in finite time except if dw/dp diverges (or p, but that would lead to a singular 
cosmology), see Ref. |208j . The physical reason is that we demand p to be a unique 
function of p, but at w = — 1 we find that p = 0. If the fluid crosses w = —1 the 
energy density p will first decrease and then increase again, while the pressure p will 
monotonically decrease (at least near the crossing). Fig. (fTO|) . It is therefore impossible 
to maintain a one-to-one relationship between p and p (our starting point), see Fig. (jlip 
(notice that the maximum of p and the minimum of p do not coincide). 

Let us have a closer look at the perturbations. The second term in the expansion 
Eq. (|182p can be re-written as 



dp 
dp 



p w _ 2 



p m{i 



^ < (184) 



where we used the equation of state and the conservation equation for the dark energy 
density in the background. We notice that the adiabatic sound speed will necessarily 
diverge for any fluid where w crosses —1. 
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Figure 10: Energy density (black solid line) and pressure (red dashed line) as functions of the 
scale factor, w crosses —1 at — 0.5, so that p is minimal at this point while p decreases 
monotonically there (but has a maximum earlier). 
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Figure 11: The pressure as function of the energy density. The graph oi p{p) shows that p is 
not a single valued function of the energy density, and that there is a point with an infinite 
slope (corresponding to the minimum of p and the divergence of c^). The point where p(p) 
has a zero slope corresponds to the maximum of p and a vanishing c^. 
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For a perfect barotropic fluid the adiabatic sound speed turns out to be tlie physical 
propagation speed of perturbations. It should therefore never be larger than the speed 
of light, otherwise the theory becomes acausal, see Refs. [2091 I210j . The condition 
Ca < 1 for (1 + w) > (our point of departure) is equivalent to 



Therefore, no barotropic fluid can pass through w = ±1 from above without violating 
causality. Even worse, the pressure perturbation defined in Eq. (|179p will necessarily 
diverge if w crosses —1 and dp ^ 0. 

• Non-adiabatic fluids 

The discussion of the barotropic fluid shows that we have to violate the constraint that 
p be a function of p alone. At the level of first order perturbation theory, this amounts 
to changing the prescription for dp which now becomes an arbitrary function of k and t. 
This problem is conceptually similar to choosing the background pressure p(t) , where 
the conventional solution is to compare the pressure with the energy density by setting 
p(t) = w{t)p(t). In this way we avoid having to deal with a dimensionfuU quantity and 
can instead set w, which has no units (up to a factor of the speed of light) and so is 
generically of order unity and often has a simple form. 

It certainly makes sense to try a similar approach for the pressure perturbation. How- 
ever, there are two relevant variables that we could compare 6p to, the fluid velocity V 
and the perturbation in the energy density Sp. Clearly it would be counterproductive 
to replace a single free function by two free functions, and it would lead to degeneracies 
between the two. Another problem is that the perturbation variables depend on the 
gauge choice. But in this case the two problems cancel each other, leading to a simpler 
solution: going to the rest-frame of the fluid both fixes the gauge and renders the fluid 
velocity physically irrelevant, so that we can now write, see Ref. jl95j : 



where a hat denotes quantities in the rest-frame. The physical interpretation is that 
cl{k, t) is the speed with which fluctuations in the fluid propagate, i.e. the sound speed. 

The rest-frame is chosen so that the energy-momentum tensor looks diagonal to an 
observer in this frame. 

Let us calculate the pressure perturbation defined by Eq. (jl86p in the conformal New- 
tonian frame, following Ref. [195J: breaking the single link between p and p amounts 
to the introduction of entropy perturbations. A gauge invariant entropy perturbation 



— <3{l + w){l~w). 



(185) 



Sp = clSp, 



(186) 



variable was found to be P ^ ^ 



w p 



, see Ref. [111 [THS]. By using 




(187) 



and the expression for the gauge transformation of 6p, Ref. [216) . 




(188) 



we find that the pressure perturbation is given by 




(189) 



^^We assume that there are no coupUngs beyond gravity between the fluids. 
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As — oo at the crossing, it is impossible that all other variables stay finite except if 
V ^ fast enough; and this is not in general the case, except if — >■ or w' — >• at 
crossing. 

To this end, we consider again the structure of the perturbation equations near it; = — 1. 
Inserting the expression for 5p into our system of perturbation equations we find 

S' = 3il + w)i,'-(^-^+3Hicl~cl)^^-^{cl~w)6 (190) 

The spoiler here is the continued presence of which we know to diverge at crossing. 

However, by construction the pressure perturbation looks perfectly fine for precisely 
the observer in the rest-frame, as Sp = c^6p does not diverge. Our prescription for the 
pressure perturbations has singled out the one frame which we cannot use for fluids 
crossing the phantom divide. The reason is that the gauge transformation relating the 
pressure perturbations in the different gauges is: 

Sp = Sp + 3Hpcl^. (192) 

If V does not vanish fast enough at the crossing then the pressure perturbation has to 
diverge in at least one frame. 

Clearly, the problem is to find a way of characterizing the pressure perturbations in a 
physical way, in Ref. |213j is shown a clear example, called the Quintom model where it 
was found that the additional contributions to the pressure perturbation diverge but in 
such a way that we end up with a finite result for Sp. In other words, even though the 
propagation speed of sound waves remains finite and constant, the additional internal 
and relative pressure perturbations lead to an apparent sound speed which diverges. It 
is only the sum of all contributions to dp which remains finite. 

At the end, models with purely adiabatic perturbations cannot cross w = — 1 without 
violating important physical constraints (like causality or smallness of the perturbations), 
it is possible to rectify the situation by allowing for non-adiabatic sources of pressure per- 
turbations. However, the parameterisation of dp in terms of the rest-frame perturbations 
of the energy density cannot be used as this frame becomes unphysical at w = —1. By 
parameterising Sp instead in any other frame the divergencies are avoided, see Ref. jl94j . 

The dark energy sector is still not well defined and several problems arise also at the 
perturbative level. The main problem here is that as far as we still believe the dark energy 
to be responsible for the acceleration we cannot demand their perturbations to be zero. 

8.3 Equivalence between general fluid and scalar field 

Here we show that at the level of first-order perturbation theory a scalar field behaves just 
like a non-adiabatic fluid with = 1. To this end we decompose the scalar fleld into a 
homogeneous mode (/){t) and a perturbation S(j)(k,t). At the background level we find 

P = ^'^' + ^('^) (193) 
p = ^c/.2-y(0) (194) 



42 



and the equation of conservation is just the equation of motion, 



2-H(/) + a^^ = 0. (195) 



The adiabatic sound speed is defined as: 

(and we remind the reader that a/a = % = aH). The perturbed energy momentum tensor 
is: 

-5T^ = 5p=\^5<f>-\^^-^ + ^5cj, (197) 

aq) 

5TI = ^p=l0<50-l02^-^J0 (198) 

do 



^ikST^ = iM7^° = -^(I^H = pV (199) 



where we wrote ^> for the gravitational potential -0 in order to avoid confusions with the 
scalar field variables (only in this subsection). 

In order to derive the rest frame sound speed of the scalar field we use equation (|189p 

Sp = ciSp+^-^{cl-cl)-pV (200) 
and express everything in terms of scalar field quantities. We find 

aq) 

-2 ( 1 IrU 1 12. dV 3aH 
a(p 

iaH ( 2a 1 dV\ ■ , , 

which after some algebraic manipulations turns into 

1 ; 1 12-^- dV ^ ^ . 3aH 



-.2 



= £2 f 1<^<50 - 1^2^ + ^Scf. + ^^sA . (202) 
\a d0 / 

Therefore — 1. 

Now let us derive the equation of motion for the scalar field perturbations from the 
perturbation equation for a perfect fiuid, Eq. (|175p . 

5=-{l + w)(e- 3*) - 3aH - w6^ (203) 

which can be rewritten as 

Sp + 3aH{p + p)S^3{/5 + p)4'~pV. (204) 
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Expressing the time derivative Sp in terms of scalar field quantities, 



Sp 



2aH ■ ■ , 



f 

+ - 
6aH 



2aH ■■, dV ■ 



(205) 



and doing likewise with the other terms, 

iaH {Sp + Sp) 
3(P + P)* 
pV 



+ 2<t)—^ + —S4> + -—r 



6aH ■ ^- , 6aH ■ . 



n2 



we can insert all these expressions into Eq. (|204p and obtain finally 



2aHS(j) + 



fc2 



dV 



(206) 
(207) 
(208) 

(209) 



which is indeed the equation of motion for 50 in the conformal Newtonian gauge (see e.g. 
Ref. [IMI)- 

8.4 Comparing different models 

We are now in the position to compare different dark energy models also at perturbation 
level. In the previous section we assumed that dark energy has no anisotropic stress which 
implies that the metric perturbation variables are the same (p — ip, see Eq. (|I70p . In general 
this is not true or better, there is no evidence to set the anisotropic stress equal to zero. This 
new quantity increases the degrees of freedom of dark energy since it depends on the intrinsic 
characteristics of the fluid itself and needs to be measured. It was shown in Ref. [214J that 
the anisotropic stress is a crucial quantity for the growth of perturbations moreover if we 
want to distinguish between dark energy and alternative gravity theories. 

We show here how the dark energy perturbations influence the dark matter and the 
metric perturbations providing an explicit example in which a general dark energy model 
can reproduce the metric perturbation of the DGP scenario; as a consequence the growth 
rate of structures is not sufficient to distinguish between dark energy and modifications of 
gravity. 

As mentioned before, the physical properties of the fluid are given by the anisotropic 
stress a and the pressure perturbation Sp (in general both can be functions of time as well 
as of k). For the dark energy all these quantities are a priori unknown functions and have 
to be measured. For the special case of dark energy due to a minimally coupled scalar field 
we have a variable w (corresponding to the choice of the scalar field potential, and fixed by 
the expansion history of the Universe), de ~ ^ ^'^"^ cr = 0, see e.g. Ref. [1961 1213] . 

The perturbations in different fluids are linked via the perturbations in the metric (p and 
■0. Introducing the comoving density perturbation A = S + iHaV/k'^, Eqs. pl7)) and ^70)) 
in standard cosmology is given by 



i 



(210) 
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Figure 12: This figure shows how the growth of tire matter perturbations depends on the 
clustering properties of dark energy. From the top downward the sound speed is = 
—2 X 10^"* (cyan dash-dotted line), = —10"'' (magenta long dashed line), = (blue 
dotted line) and = 1 (red dashed line) . For comparison we also plot the growth factor of 
the DGP model (black solid line). 

where the sum runs over matter and dark energy in our case. 

The quantity we are interested in is the growth factor G = Am/a which parameterises 
the growth of structure in dark matter. The growth factor is normalized so that G = 1 for 
a <C 1 (using that cx a during matter domination and on sub- horizon scales). We assume 
that G is an observable quantity (even though of course large scale structure surveys observe 
luminous baryonic matter, not dark matter, adding yet another layer of complications). 

The growth factor is not uniquely determined by the expansion history of the Universe 
(and hence w of dark energy). Although the main effect of the dark energy is to change H , 
leading to G < 1 at late times, there is an additional link through the gravitational potential 
■0. Different dark energy perturbations will lead to a different evolution of ip, which can 
modify the behavior of G. Conventionally dark energy perturbations are neglected because 
are considered to be unimportant, e.g. Ref. [219|. This is a good assumption for scalar field 
dark energy where the high sound speed prevents clustering on basically all scales. However, 
a small sound speed de ~ is not excluded. Indeed, it could also vary in time. We show in 
Fig. p2p how the growth factor of the dark matter changes in response to large dark energy 
perturbations 0. 

What happens is that, as we decrease the sound speed, dark energy is able to cluster more 
and more. The increased dark energy perturbations lead to enhanced metric perturbations. 
The dark matter in turn falls into the potential wells created by dark energy, leading to an 
increase of the growth factor. Although clearly G is not uniquely determined by wde, we 
notice that it always increases as we decrease (at least as long as the linearized theory 

is applicable). Looking at the evolution equations (|180p and ()18ip for cr = (<J=> cj) — ip) 
we see that the response of the fluids to the metric perturbations is governed by the sign of 

^®We emphasize that we discuss only how dark energy perturbations can modify the behavior of dark 
matter, without taking into account limits from observations. 
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1 + w. Non-phantom dark energy clusters therefore in fundamentally the same way as the 
dark matter and can only increase the growth of matter relative to the case of negligible 
dark energy perturbations. 

So although dark energy perturbations can influence the growth factor of dark matter, 
they only seem capable of enhancing it. In Fig. (|12p is also shown the prediction for the 
growth factor in the DGP model from Ref. [217] . and it is smaller than the one of a smooth 
dark energy component. Therefore, dark matter clusters weaker in DGP model than in any 
other model considered so far. Clearly, a sensitive experiment will be able to distinguish the 
different curves and rule out the wrong model. 

A natural question is: is it possible to have a dark energy model able to mimic DGP 
model also at perturbation level? 

To this end we need to take a closer look at the DGP model. 

An important aspect of DGP and other brane- world models is that the dark matter does 
not see the higher-dimensional aspects of the theory as it is bound to the three-dimensional 
brane. Its evolution is then the same as in the standard model. The modifications appear 
only in the gravitational sector, represented by the metric perturbations. 

The metric perturbation in DGP can be written as, see Refs. [2171 1218| 



k'cj) = -47rGa^ ( 1- — )p„A™ (212) 
fcV = -47rGaMl + -!-)p™A™ (213) 



3/3 



where the parameter (3 is defined as: 



13 = 1- 2r,H 1^1 + ^ j = 1 -I- 2r,HwuE. (214) 

Dark matter does not care if the metric perturbations are generated (in addition to its own 
contribution) by a modification of gravity or by an additional dark energy fluid. Its response 
to them is identical. Or to put it differently, if dark energy and dark matter together can 
create the ip and ip of Eqs. (|212p and pi3p then the growth factor (and indeed all other 
cosmological observables) will be the same as in the DGP scenario. 

We see immediately that in order to generate these metric perturbations we will need 
to introduce an anisotropic stress since (p ^ tp. This seems to be a very generic property 
of modified gravity that is also present in f{R) models, see Ref. [220' and has been noticed 
before. We plot in Fig. p3p again the growth factor for scalar field dark energy and the DGP 
model, but now also a family of dark energy models with non-vanishing anisotropic stress 
a. We notice that these models can easily suppress the growth of perturbations in the dark 
matter for cr < and mimic the behavior of the DGP model. 

Formally we can recover the DGP metric perturbations by choosing 

'^DE = T— A„, (215) 

9/3(1 -I- wde) Pde 

for the anisotropic stress of the dark energy, if we can also generate dark energy perturbations 
with ^ 

PDeAdE = -—Pm^rn- (216) 

We notice that these are very large dark energy perturbations. Indeed, if we keep = 1 and 
set cr to the expression piSp we suppress the growth of the matter perturbations too much, 
see Fig. ([T^ . Since /3 < the large dark energy perturbations of Eq. (|216p then increase the 
matter clustering back to the DGP value. 
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Figure 13: In this figure we show how the anisotropic stress of the dark energy afi^ects the 
growth of the dark matter perturbations. The red dashed fine corresponds to scalar field 
dark energy with — 1 and a — 0. The dotted blue line shows how the dark matter 
growth factor decreases for a constant ctde = —0.1. The long-dashed magenta line uses the 
theoretical anisotropic stress of Eq. (|215p with = 1, which suppresses the growth of the 
matter perturbations too much. Finally, the dash-dotted cyan line (nearly on top of black 
solid DGP line) uses the same ctde but sets the pressure perturbation of dark energy to 
6p = (1 + w)pa in its rest frame. 
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The required size of the dark energy perturbations in itself is no problem, as we can lower 
the sound speed and even make it negative. However, while for a = we were not able to 
decrease Am with the help of dark energy perturbations, we find that with a large, negative 
anisotropic stress we are unable to increase it. The required anisotropic stress is far larger 
than the gravitational potential ^p, and it starts to be the dominant source of dark energy 
clustering in Eq. (|176p . As it enters with the opposite sign it now leads to anti-clustering of 
the dark energy with respect to dark matter which feels only ip (i.e. dark matter overdensities 
are dark energy voids). There is still enough freedom in the choice of a to match the growth 
factor very precisely, but if we could measure both (f> and separately then we could detect 
the differences between the two models. 

Is it really not possible to match both ip and of the DGP model within a generalized 
fluid dark energy model? Yes, it is: The metric perturbations have two degrees of freedom, 
and we do have two degrees of freedom of the dark energy to adjust, a and 5p. Allowing 
free use of the pressure perturbations, we can choose them for example to cancel the direct 
effect of a onto the dark energy perturbation in Eq. (|176p by setting 6p = (1 + w)pa. 
This reverses the sign of Ade, and minor adjustments to the pressure perturbations can 
then provide the required match to A,„. For the cyan dash-dotted curve in Fig. [T3] we set 
5p = {1 + w)p(j + 3Hac^pV/k^, i.e. we cancelled the contribution of the anisotropic stress in 
the dark energy rest frame. This provides a very good solution to Eqs. (|215p and (|216p during 
matter domination. It is easy to improve the solution to the point where it is impossible 
to distinguish observationally between the DGP scenario and this generalised dark energy 
model. 

We can be even more general since the reader might have thought that the previous 
discussion worked only for the DGP model. 

Let us assume that the (dark) matter is three-dimensional and conserved, and that it does 
not have any direct interactions beyond gravity. We assume further that it and the photons 
move on geodesies of the same (possibly effective) 3+1 dimensional space-time metric. In 
this case we can write the modified Einstein equations as 

X^, = -SttGT^, (217) 

where the matter energy momentum tensor still obeys T^.^, — 0. While in GR this is a 
consequence of the Bianchi identity, this is now no longer the case and so this is an additional 
condition on the behavior of the matter. 

In this case, we can construct Y^i, = X^i/ — G^,/, so that G^ti/ is the Einstein tensor of 
the 3+1 dimensional space-time metric and we have that 

Gpi/ — — SttGT^i/ — Yf^jy. (218) 

Up to the prefactor we can consider Y to be the energy momentum tensor of a dark energy 
component. This component is also covariantly conserved since T is and since G obeys the 
Bianchi identities. The equations governing the matter are going to be exactly the same, by 
construction, so that the effective dark energy described by Y mimics the modified gravity 
model. 

By looking at Y we can then for example extract an effective anisotropic stress and an 
effective pressure perturbation and build a dark energy model which mimics the modified 
gravity model and leads to exactly the same observational properties, see Refs. |214j and 

m- 

This is both good and bad. It is bad since cosmology cannot directly distinguish dark 
energy from modified gravity )although there could be clear hints, e.g. a large anisotropic 
stress would favor modified gravity since in these models it occurs generically while scalar 
fields have tt^ = 0, see for instance Ref. |221| ). However, it is good since there is a clear target 
for future experiments: Their job is to measure the two additional functions describing Y as 
precisely as possible. 
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We will show in the next sections how baryon acoustic perturbations and weak lensing 
can help to constraint H, (j) and ip. 

9 A possible test 

We have shown in the previous section that the growth factor is not sufficient to distinguish 
between modified gravity and generalised dark energy, even if the expansion history (and so 
the effective equation of state of the dark energy) has been fixed by observations. Dark energy 
models can match the metric perturbations completely so that cosmological observations 
cannot distinguish between the two possibilities, as we have shown explicitly for the DGP 
case. 

Although the construction of a matching dark energy model for the DGP case may seem 
very fine tuned (and it is), we were more concerned with the question to what degree this 
is possible at all. Just measuring a growth factor that does not agree with scalar field dark 
energy is not sufficient to rule out dark energy and General Relativity. But clearly, if the 
expansion history and the growth of matter perturbations were to match those predicted 
from a physically motivated and self-consistent modified gravity model, a statistical analysis 
would rule out a fine tuned dark energy model. 

The big question now is: is there something that we can do? Among all the models we 
have at the moment to describe the acceleration of the Universe, can we rule out some of them 
from theoretical and/or observational considerations? The answer is somewhat ambiguous. 
All the dark energy models that we can found in the literature are completely arbitrary (also 
those we reported in this review which are considered for faith as the most accepted ones) 
until fundamental physics is able to justify the form of the potential^. 

As we saw, the ACDM is the simplest model to explain the acceleration and it is seems to 
be also the most accepted one because it is confirmed by all the observations we now possess. 
However, as we have also shown previously, it is always possible to have a scalar field which 
manifest an equation of state parameter close to —1 and we will never know the difference. 

However, this should not discourage the reader because we can always require some 
properties that our model has to respect!^. We can consider it as a test that a particular 
model needs to pass in order to be considered as a good candidate. We should always have 
in mind that our knowledge on dark energy is still poor and all the points listed below can 
be still debatable. 

• If we are dealing with a real form of matter /energy 

We list here some properties that we think are fundamental for our model to be accepted. 
Should our theory follow: 

1. Lorentz invariance?: this means that the theory has to be not only covariant under 
transformations but also that the only absolute quantities possible in the theory are 
pure constants and they are the same in all the Universe. 

2. Energy density?: as mentioned in the previous sections, the energy density of a 
dark energy candidate should be in agreement with observations and the coincidence 
problem must be solved. In other words we should be able to solve the fine tuning 
problem. 

3. Close to cosmological constant?: at the moment all the data we possess favour dark 
energy models whose equation of state parameters are close to —1. Probably we should 

^^If we are dealing with a scalar field being it associated to quintessence, K-essence etc... 
However these are still open questions. 
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think of it as a clear sign that regardless of what is the responsible for acceleration, a 
model should manifest a. w ^ —1 (effective or not). However, it is still not clear from 
theoretical basis if the equation of state parameter can be less than —1 (i.e. phantom). 
From the cosmological point of view there is no problem, however classical general 
relativity imposes energy conditions restricting the range of available values of < 1, 
see Ref. pii]. 

4. Ghost?: when a theory has a ghost it means that the kinetic term in the Lagrangian 
has a wrong sign; the main problem comes when we want to quantize the theory as it 
is not possible to define a ground state. So, in order for a theory to be considered as a 
good theory should be ghost-free. 

5. Causality?: each theory which is also Lorentz invariant should not send information 
with a speed larger than the speed of light: c = 1. In practice, we can say that the 
sound speed (speed with which perturbations, and hence information, propagate in the 
fluid) should not be bigger than unity: < 1 (and possibly also positive in order to 
prevent perturbations to blow up). 

6. Evolution of density perturbations?: a dark energy model should lead to a proper 
evolution of the mass density field consistent with the observations: CMB, matter 
power spectrum, peculiar velocity and weak lensing power spectrum. 

• If we are dealing with a modification of gravity 

An obvious requirement is that the theory should be well defined and it has to follow 
the same criteria listed above; moreover, should our model also follow: 

7. Proper evolution of different epochs?: a modified gravity model should account 
for a correct evolution of the scale factor in the matter dominated era (as required from 
structure formations) . 

8. Local gravity test?: we are in the position to accept that General Relativity works 
extremely well in our local Universe, hence a particular model of modified gravity 
should pass the solar system test, i.e. / (R) R — 2A. 

According to the requirements listed above we can examine the models we presented in 
this review and see whether or not they pass the tests; in Tab. 1 are presented a subset 
of dark energy candidate, when the candidate satisfies the corresponding requirement the 
symbol /is used whereas the symbol x is used instead. If the symbol /is followed by the 
symbol ! this means that the model partially satisfies the criteria (or it is still not clear 
whether or not it is applicable). Clearly a candidate can be ruled out if contains at least one 
X. 
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Candidate 


tests 




1 


2 


3 


4 


5 


6 


7 


8 


A 




X 




/ 


/ 






•J 


Quintessence 


/ 


X 


/ 


/ 


/ 


/ 


/ 


/ 


k-essence 


/ 


/! 


/ 


/ 


X 


/ 


/ 


/ 


Phantom 


/ 


X 


/ 


X 


/ 


/ 


/ 


/ 


Chaplygin 


/ 


/! 


/ 


/ 


X 


/! 


/ 


/ 


Coupled DM-DE 


/ 


/! 


/ 


/ 


/ 


/! 


/ 


/ 


f{R) 


/ 




/! 


/ 


/ 


/! 


/ 


X 


DGP 


/ 




X 


X 


X 


/! 


/ 


X 



Table 1: The symbol /is used when the models satisfies the corresponding requirement; the 
symbol x is used when it does not satisfy the criteria. The symbol /! means that the model 
partially satisfies the criteria (or it is still not clear whether or not it is applicable). The 
empty boxes mean that it is not required to the model to satisfy the corresponding criteria. 

10 The importance of being phenomenological 

In this section we introduce again some concepts that we already treated before and we 
apologize if sometimes we look repetitive but we believe that recalling concepts may set the 
discussion clearer. 

An alternative approach to understand the observed late-time accelerated expansion of 
the Universe is to construct general parameterisations of dark energy, in the hope that 
measuring these parameters will give us some insight into the mechanism underlying the 
dark energy phenomenon. 

A successful example for such a phenomenological parameterisation in the dark energy 
context is the equation of state parameter of the dark energy component, w = p/ p. If we 
can consider the Universe as evolving like a homogeneous and isotropic FLRW Universe, 
and if dark energy is not coupled to anything except through gravity, then w{z) completely 
specifies its evolution: as it was mentioned before, dark energy (or anything else) is described 
by the homogeneous energy density pde and the isotropic pressure pde, corresponding to 
the Tg and Tl elements respectively in the energy momentum tensor in the rest-frame of 
the dark energy. Any other non-zero components would require us to go beyond the FLRW 
description of the Universe. 

As it was mentioned before, the equation of state parameter can be measured from SNIa. 
However, to include other observables like CMB, Weak Lensing and Matter power spectrum 
(we will come back on these topics) we need to improve our description of the Universe, by 
using perturbation theory. If we work in the Newtonian gauge, we add two gravitational 
potentials </) and ip to the metric, see Section [8T2l They can be considered as being similar to 
H in that they enter the description of the space-time (but they are functions of scale as well 
as time). Also the energy momentum tensor becomes more general, and p is complemented 
by perturbations 5p as well as a velocity Vi . The pressure p now can also have perturbations 
6p and there can further be an anisotropic stress a. 

The reason why we grouped the new parameters in this way is to emphasize their role: 
at the background level, the evolution of the Universe is described by H, which is linked to 
p by the Einstein equations, and p controls the evolution of p but is a priori a free quantity 
describing the physical properties of the fluid. Now in addition there are and ip describing 
the Universe, and they are linked to 6pi and Vi of the fluids through the Einstein equations. 
6pi and a in turn describe the fluids. Actually, there is a simplification: the total anisotropic 
stress a directly controls the difference between the potentials, (f> — ip- 

The gravitational potential is given by Eq. piOp . All the fluids with non-zero perturba- 
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tions will contribute to it. If we assume for instance that dark energy does exist we cannot 
demand that its perturbations are too small. In addition it might be that we are dealing with 
some modified gravity model or something more complicated, see for instance Eq. (|212p : in 
order to account for all these effects we parameterise the Poisson's equation as: 

= -AnGa^Q {k, a) p™A„, (219) 

where A = 5 + 3aHV/k'^ is always the comoving density contrast. 

Here Q{k,a) is a phenomenological quantity that, in General Relativity, is due to the 
contributions of the non- matter fluids (and in this case depends on their 5p and cr). But it 
is more general, as it can describe a change of the gravitational constant G due to a mod- 
ification of gravity. It could even be apparent: If there is non-clustering early quintessence 
contributing to the expansion rate after last scattering then we added its contribution to the 
total energy density during that period wrongly to the dark matter, through the definition 
of Om- In this case we will observe less clustering than expected, and we need to be able to 
model this aspect. This is the role of Q{k, a). 

For a dark energy model we have to admit an arbitrary anisotropic stress a (see Eq. (|211l) ) 
and we use it to parameterise -0: 

i/; =[l + r]{k,a)](t). (220) 

There is no sign for a non-vanishing anisotropic stress beyond that generated by the free- 
streaming of photons and neutrinos. However, it is expected to be non-zero in the case of 
topological defects, see Ref. |227j or very generically it is expected to be non-zero for modified 
gravity models (see Eq. (|213p V With Q and rj (usually called in the literature gravitational 
slip) then we can recover all the possible models. 

However, it is not always possible to have an analytic expression for the parameters above 
mostly because the equation of perturbations are too complicated to be solved. Of course, 
it is always possible to solve them numerically. However, analytical results allow a much 
better insight and also an easy way to see how the behavior of some observables changes 
as a function of the parameters. For instance, for all the class of dark energy models the 
parameter Q (fc, t) can be thought as: 

Q(fc,i) = l + ^^£^. (221) 
111 Ref. |228j can be found an analytical expression for Q {k,t): 

" ' ^"(l + H . " ...... (222) 



here ri{k,t) = because it was assumed zero anisotropic stress for dark energy. This ex- 
pression is evaluated under the assumption that both and w are constant. This latter 
assumption is usually violated, but the results should nonetheless allow insight into the be- 
havior of the perturbations. For models where the quantities vary only slowly with time, it 
is expected that the results will still hold in an averaged sense, due to the indirect nature 
of most observations. The big advantage of making these assumptions is that it allows us 
to solve the perturbations analytically under the additional condition of matter domination, 
leading to surprisingly simple results. In addition, the resulting expression for Q is surpris- 
ingly accurate even for late times. The reason is that both fluids, dark energy and dark 
matter, respond similarly to the change in the expansion rate so that most of the deviations 
cancel, see Eq. p21|) . 

We would like now to connect the parameterisation introduced before to the observables. 
A very promising experiment regards weak lensing (we will come back later in more details), 
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the distortion of foreground galaxies will help us to constrain cosmological parameters, like 
dark energy parameter of equation of state, with extraordinary accuracy. 

The weak lensing effect is proportional to the lensing potential $ = (p + ip, it is then 
useful to consider a new parameterisation: 

= Q(2 + r;)5^|^A™ = 2S^^|^A„,. (223) 
2a 2a 

Correspondingly, the total lensing effect is obtained by multiplying the usual equations by a 
factor S = (1 + ri/2)Q. 

In linear perturbation theory all k modes evolve independently, so that we can decompose 
the dark matter density contrast as: 

A„, = aG (a, k) A„ (fc) . (224) 

Here A^ (fc) = A„i (a = 1, fc) determines the matter power spectrum today, P{k) = | Am(fc)P- 
G (a, fc) is the growth factor which is scale-independent in the usual ACDM cosmology. How- 
ever, the contribution from the dark energy perturbations induces a scale dependence, as we 
showed in the previous section. 

However, this is not the only way the evolution of dark matter perturbations is affected 
by the evolution of dark energy perturbations. To see this we can write the weak lensing 
potential as: 

fc2$ = _3iJ (af a^E (a, k) (a) G (a, fc) A,„ (fc) (225) 

where we used Eq. (|224|) . 

Hence, the lensing potential contains three conceptually different contributions from the 
dark energy perturbations: 

• The direct contribution of the perturbations to the gravitational potential through the 
factor S. 

• The impact of the dark energy perturbations on the growth rate of the dark matter 
perturbations, affecting the time dependence of A,„, through G (a, fc). 



A change in the shape of the matter power spectrum P{k), corresponding to the dark 
energy induced fc dependence of A^. 



In the standard ACDM model of cosmology, the dark matter perturbations on sub- horizon 
scales grow linearly with the scale factor a during matter domination. During radiation 
domination they grow logarithmically, and also at late times, when the dark energy begins 
to dominate, their growth is slowed. The growth factor G = /^m/a is therefore expected to 
be constant at early times (but after matter-radiation equality) and to decrease at late times. 
In addition to this effect which is due to the expansion rate of the Universe, there is also the 
possibility that fluctuations in the dark energy can change the gravitational potentials and 
so affect the dark matter clustering. 

In ACDM G can be approximated very well through 



G{a) = cxp ^^J d\na( 



{n„,{ay'-l)\ (226) 



where 7 « 0.545. There are two ways that the growth rate can be changed with respect to 
ACDM: Firstly, a general w of the dark energy will lead to a different expansion rate, and 
so to a different Hubble drag. Secondly, if in the Poisson equation pi9p we have Q ^ I then 
this will also affect the growth rate of the dark matter, as will r/ 7^ 0. We therefore expect 
that 7 is a function of w, rj and Q. 
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For standard quintessence models we have 77 = and Q w 1 on small scales as the scalar 
field does not cluster due to a sound speed — 1. In this case 7 is only a weak function of 
w. A constant 7 turns out to be an excellent approximation for some coupled dark energy 
models as well, see Ref. [239] . and for some modified gravity models, see Ref. |231j . A 
similar change in 7 can also be obtained in models where the effective dark energy clusters, 
see Refs. pTllMI] . 

However, we would like of course to be more general; on sub- horizon scales it can be 
assumed that k'^ip dominates over (f)' as the source of matter clustering in the equation of 
perturbations, also because the potentials are normally slowly varying. The term with ip is 
now multiplied by {I + r])Q and following the discussion in Ref. j241| we have: 

■^-'^^^ (227) 

1 - ilm[a) 

and the asymptotic growth index becomes: 

3{l-w^-A{Q,f^)) 

^ — 5^6^^;;^ — ■ (22^) 

Although these relations are useful to gain physical insight and for a rough idea of what 
to expect, for an actual data analysis one would specify Q and 77 and then integrate the 
perturbation equations. 

Then, the set of parameters {7, S} is completely general and it can be used also for 
modified gravity models as we can reconstruct a set of effective parameters also for this class 
of models. 

Therefore, there is a clear target for future experiments: together with w, only two extra 
quantities {Sp and a or equivalently 7 and S) are needed to be measured in order to span 
the complete model space for both modified gravity and dark energy models. We can find 
first observational tests and future constraints on the sound speed and the anisotropic stress 
in Refs. [242], [243], [244], [24g and [246] . 

11 The new observational tests 
11.1 Baryon acoustic oscillations: BAO 

Let us now consider in more details the standard rulers introduced in the section 15.41 If we 
observe a transverse comoving scale Ai which subtends an angle 6, then the angular diameter 
distance is 

Di Xi/ {1 + z)e (229) 

where the subscript indicates a given cosmology. In a different cosmology, the angular diam- 
eter distance will be given by 

D2 = X2/il + z)9, (230) 

in brief we can say that the scale has to change in order to keep the same subtended angle 
at the same redshift. Comparing Eq. p29|) and Eq. ()230|) it can be seen that the ration D/X 
is constant for any given angle. We can write everything in terms of wavenumber; therefore 
if we refer with the subscript r to the reference cosmology, then we have that for any other 
cosmology the transverse wavenumber is given by: 

fc_L = fcr,±^- (231) 

The same discussion can be done with the comoving scale extending a long the line of sight 
from two different redshifts, say zi and Z2- The scale is given by: 

dz 

'=Wz) (2^2) 
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where dz = Z2 — zi. In order for this scale to be seen at the same dz the product XH{z) has 
to remain constant when changing the cosmology. As before, we can have a relation for the 
radial mode: 

.„=^„|^. (233) 

The relations above can be applied to any mode. Moreover, every mode can be written in 
term of a reference mode kr with a dependence on the cosmological parameters: 

k= (kl + kl^ ^'"^ = Rkr (234) 

where 

^ - h;d ^^^^^ 

and yit = k • x/A: is the direction cosine along the line of sight, here x is the unit vector parallel 
to the line of sight. 

We showed that the wavenumber changes with the cosmology, this implies that if the 
matter power spectrum is isotropic in the reference cosmology it will be anisotropic for any 
other cosmology since k\\ and k± scale differently: this is called the Alcock-Paczynski effect, 
seeRef. [233] . 

Moreover, the power spectrum is proportional to the volume V in which we measure the 
perturbations, hence if cosmology changes then we need to take into account also the change 
of the volume. Let us assume that we measure the power spectrum within a solid angle 0^ 
and a shell of width dz, then the volume is: 

V = 9\^dr = e^'^dz=^^-^dz (236) 
H{z) H{z) 

where dr = dz/H{z) and D'^{z) — 9'^r'^{z). Following the same argument used before, we 
can express the volume in terms of the reference volume, it reads: 

The matter power spectrum is simply a product of the spectrum of primordial fluctuations 
and the modifications of those fluctuations in the later epoch. Linear perturbation theory 
fixes the matter power spectrum in comoving coordinates and changes only the amplitude 
as the structures evolve (once matter-radiation equality is fixed). The growth factor G (z) 
rescales the amplitude of the fixed matter power spectrum to account for the growth of 
structures from the recombination to a redshift z. Then the matter power spectrum can be 
further simplified as: 

P{z,k) ^G^{z)P{k) (238) 

where G{z) = 5m{z) / SmiO) is called growth factoi0, where S is the energy density contrast 
of the matter field. 

However, the observable in the Universe is the galaxy over-density which is assumed to 
trace the underlying matter distribution through a function called bias factor. In principle 
this quantity could be arbitrary; it could depend either on scale and time. Usually it is 
assumed that the bias on large scale is independent on scale hence, in the matter power 
spectrum, this term appears as a multiplicative factor which modulate the overall amplitude 
of the galaxy power spectrum, b'^{z). It is worth noticing that there are works which pointed 
out that this assumption is too strong and it may be wrong, see Ref. |234| . 



^^This is certainly true if dark energy is accounted for only at background level. The growth factor can 
also be a function of space if dark energy can cluster as pointed out before. 
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Then, an observer can only measure the galaxy power spectrum in redshift space, which is 
distorted compared to the power spectrum in the real space, where the theoretical predictions 
were made. Redshift distortions are angle dependent distortion on the power spectrum caused 
by the peculiar velocities of galaxies. In linear theory these will affect structures along the 
line of sight; i.e. for a structure which is isotropic in the real space, an observer will measure 
more power in the radial direction than in the transverse direction. On large scales, these 
distortions follow a simple form and it appears as a multiplicative angle dependent factor in 
the power spectrum, (l + , see Ref. |235j . The faction (3 is defined hy f5 — f /h where 

/ is the growth rate of matter perturbations defined as: 

Putting all the effects together we finally have the observed galaxy power spectrum: 

Pobs{z. K) = ^l^§^GHz)bHz) (1 + /3m')' Por{k) + Pshotiz) (240) 

where we added the shot noise Ps{z) which depends only on redshift and it is a Poissonian-like 
noise coming directly from the number of galaxy counted in the survey volume. 

Another layer of complication is the non linear evolution of the matter density fields. As 
time goes by matter density fields start to interact with each other breaking the linearity, 
linear theory does not describe correctly the evolution of structures if observations are made 
in the local Universe. As a consequence, this will increase the power spectrum at small 
scales (large wavenumber) starting to become an important term; non linear effects will 
automatically smear out the acoustic peaks and all the information is lost. Therefore, in 
order to observe more peaks, telescopes need to measure galaxies power spectrum at high 
redshift where the non linear effects are still small. 

What does it mean to measure the BAO scale in a survey covering different redshifts? Let 
suppose we have a survey from redshift zi up to Z2 divided in different bins of width Az. Let 
us also assume we measure the correlation function as galaxy pairs in different bins. Then, 
our estimate of the correlation function will be the average of the correlation functions in 
different bins. BAO on the power spectrum correspond to a bump in the correlation function. 

Measuring the BAO in the galaxy power spectrum will help us to say more about the dark 
energy parameters, such as: its energy density VLde and the parameter of equation of state w, 
as it will measure the Hubble parameter and angular diameter distance. Moreover, the galaxy 
power spectrum depends on the growth factor which happens to be an important quantity 
nowadays. A test of the growth factor would be important both as a consistency check for the 
standard cosmological model because G is also determined by the Hubble parameter and as 
a constraint on non-standard models like e.g. modify gravity, see section [8^ In fact, models 
that modify the Poisson equation will also generically modify the perturbation equation for 
the matter density contrast 5. 

We report here the constraints on the cosmological parameters from the latest experi- 
ments. In Fig. (fT4l) the confidence region for the parameters Vim and is plotted. It is 
important to notice that the CMB information is important to reduce the errors especially 
on the matter density parameter but the errors on the equation of state parameter wde are 
basically unaltered. 

The reason is that at moment we have survey only at low redshift [z < 0.5) whereas the 
CMB peaks give a measurement at redshift z ^ 1100 but some parameter, like wde, at that 
redshift are very diluted. 

In Fig. (|15p the constraints for flm, wde ai'e shown including the set of supernovae given 
in Ref. [211. 
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Figure 14: Likelihood surface in the (^rmWDE) plane assuming a flat space with constant 
dark energy equation of state, for spectra in the 2dF and SDSS redshift survey. The shaded 
regions show the likelihood for the redshift survey alone. The solid contours are calculated 
by modelling the CMB sound horizon scale and the dashed contours by including the CMB 
peak position measurement. From Ref. [247j 




Figure 15: As Fig. (fT4|) . but now additionally using the SNIa in the Likelihood calculation. 
The shaded region, dashed and solid contours were calculated using the BAO based measure- 
ments described in the caption to Fig. The dot-dashed contours show the likelihood 
surface calculated from just the SNLS data. 
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11.2 Weak Lensing: WL 

The idea of focusing on WL comes from two reasons: first, contrary to eg supernovae or 
baryon-oscillation tests, WL makes use of both background and hnear perturbation dynam- 
ics, allowing to break degeneracies that arise at the background level (this is particularly 
important for testing modified gravity); second, several groups are planning or proposing 
large WL experiments in the next decade (e.g. Refs. |223[ 12241 1225j ) that will reach the sen- 
sitivity to test cosmology at unprecedented depth and it is therefore important to optimize 
the science return of these proposals. We find it inappropriate to write the whole calculation 
for the lensing effect, so here we report the main features. Excellent reviews on lensing can 
be found in Refs. PSII^I^ . 

Light rays are deflected when they propagate through an inhomogeneous gravitational 
field. It was Einstein theory which elevated the deflection of light by masses from a hypothesis 
to a firm prediction. Assuming light behaves like a stream of particles, its deflection can be 
calculated within Newton's theory of gravitation, but General Relativity predicts that the 
effect is two times larger. The first experiment was in the f9f9 when the deflection of a 
light ray, coming from a star, due to the gravitational field of the Sun was measured. The 
confirmation of the larger value was the most important step towards accepting General 
Relativity as the correct theory of gravity. 

Cosmic bodies more distant, more massive, or more compact than the Sun can bend light 
rays from a single source sufficiently strongly so that we can also have multiple images Tidal 
gravitation fields lead to differential deflection of light bundles. The size and shape of the 
sources are therefore changed. Since photons are neither emitted nor absorbed in the process 
of gravitational light deflection, the intrinsic surface brightness of lensed sources remains 
unchanged. Changing the size of the cross section of a light bundle therefore changes the flux. 
Since astronomical sources like galaxies are not circular, this deformation is generally very 
difficult to identify in individual images. In some cases, however, the distortion is so strong 
that it can be recognized (this is the case of Einstein rings and arcs in galaxy clusters). Such 
strong effects are not very common in nature and most of the time we are dealing with very 
weak distortions which, in turn, are difficult to see. Although weak distortions in individual 
images cannot be detected, the net distortion averaged over an ensemble of images can still 
be detected. This is called weak gravitational lensing. 

There are two reason why gravitational lensing has become a powerful tool in modern 
cosmology: 

• The deflection angle of a light ray is determined by the gravitational fields (^/; + </)) of 
the matter distribution along its path. According to General Relativity, the gravita- 
tional fields are determined by the energy momentum tensor of matter distribution. 
So information about the matter density distribution can be obtained. Moreover, the 
light distortion is sensible to the overall mass distribution and not only to the lumi- 
nous matter, hence weak lensing is not affected by biasing uncertainty (as it was for 
the BAG). 

• Once the deflection angle is given, gravitational lensing can be easily reproduced. Since 
most of the lens system involves sources at moderate and high redshift, lensing can 
probe the geometry of the Universe. 

The two main distortion effects here that can be used in weak lensing to probe the 
statistical properties of the matter distribution between the observer and an ensemble of 
distant sources are: shear and magnification. 

The images are distorted in shape and size. The shape distortion is due to the tidal 
component of the gravitational field, described by the shear, whereas the magnification is 
given by both isotropic focusing caused by the local matter density and anisotropic focusing 
due to the shear, hence lensing changes the apparent brightness of a source. 



58 



For cosmological use, the weak lensing equation takes a simple form because of different 
assumption that can be made: 

• Density perturbations are well localized in an homogeneous and isotropic background 
(each perturbation can be surrounded by a spatially flat neighborhood). 

• The Newtonian potential of the perturbations is small and typical velocities are much 
smaller than the speed of light. 

Light rays are deflected by gravitational potential due to a mass distribution, then an 
effective surface mass density can be defined (and it will define the impact parameter). 
The power spectrum of the effective surface mass density is closely related to the power 
spectrum of the matter fluctuations (the central physical object in cosmology). Any two 
point statistics of cosmic magnification and cosmic shear can be then expressed in terms of 
the effective convergence power spectrum. 

Briefly we can say that the convergence k is related to the first derivative of the de- 
flection angle a which is proportional to the first derivative of the gravitational potential. 
The effective convergence Ke// then involves the Laplacian of the potential which is the 
Poisson's equation. The effective convergence along a light ray is therefore an integral over 
the density contrast along the light path, weighted by a combination of comoving angular 
diameter distances (sources and lens at different redshift) . Of course we are interested in the 
statistical properties of the effective convergence, especially the power spectrum. According 
to the definition above the convergence power spectrum P^ii) is a weighted integral of the 
matter power spectrum; then any weak lensing experiment will give information on the mass 
distribution (more correctly the gravitational potential) in the Universe. 

The convergence weak lensing power spectrum (which in linear regime is identical to the 
ellipticity power spectrum) is a linear function of the matter power spectrum convoluted with 
the lensing properties of space. In ACDM cosmology it can be written as, see Ref. |250| : 



where £ denotes the angular multiple, P„/[P/(fc, z)] is the non-linear matter power spectrum 
at redshift z obtained correcting the linear power spectrum Pi{k^ z) and the weight function 
Wi is the efficiency for lensing a population of source galaxy and it is determined by the 
distribution function of source and lens galaxy. 

However, there are some difficulties. Besides the problem of the linear correction, there 
are several source of systematics that can affect the outcome of the experiment. A good 
example of a work dealing with many of these details and how to control their effects is 
reported in Ref. [249]. 

The statistical uncertainty in measuring the shear power spectrum is, see Ref. |251) : 



where 7i„i is the rms intrinsic shear (arising from the intrinsic ellipticity of the galaxies and 
the measurement noise), rie// is the number of galaxy per steradians belonging to the i-th 
bin. 

WL test will be the new frontier of the experiments and powerful probes on the evolution 
of the dark energy parameter. As it can be noticed the convergence power spectrum depends 
on the matter power spectrum. As often pointed out, the growth of the matter perturbations 
depends strongly on the behavior of the dark energy and its perturbations. If we refer to 
the parameter introduced in the section [TOl we immediately realize that constraints on the 
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Figure 16: The noii linear power spectrum for ACDM and DGP model for one single bin. 
The red (central) solid line and the blue dashed line show the convergence power spectrum 
for the ACDM and the DGP model respectively. The error-bars represent the noise errors 
on the ACDM convergence power spectrum in logarithmically spaced bins. From Ref. |222) . 

set of parameters {w, 7, S} can be inferred or equivalently on the intrinsic characteristic of 
dark energy, i.e. {w,dp,a}. 

This is not the whole story, we have to remember the set of parameters {w,7,S]} is a 
generic way to express our ignorance. We should have clear in mind that these parameters 
do not necessarily refer to a real dark energy fluid but they can be also effective parameters. 
Having an good measurement of this set of parameters will enable us to say something more 
on the nature of our Universe. 

In Fig. (|16p is reported the weak lensing spectrum for ACDM with the noise due to the 
intrinsic ellipticity and for comparison the DGP spectrum for a survey such as DUNE/Euclid, 
assuming a fraction of sky fsky = 1/2 and an effective number density of Uef / = 35 galaxy 
per sq. arcmin. It can be seen that the DGP spectrum is outside the noise at low ts. 

So far, the existing cosmic shear estimation have provided interesting and complementary 
constraints on ri„j and on the spectrum amplitude (Tg but not directly on the dark energy 
parameters. These are shown in Fig. PT)) . taken from Ref. |248| . 

An increasing number of papers are giving forecasts on cosmological and on dark energy 
parameters, putting in evidence the outstanding capability of such experiments to constrain 
parameters with an extremely high precision, and possibly capable of ruling out different 
models, see Refs. [2221 ESS ElU [2541 [2551 [256] . 

In Fig. ([T8|) are shown the statistical constraints that the proposed SNAP mission could 
achieve, combining SN, WL observation and results from CMB Planck mission. 

11.3 Other probes 

Besides the methods discussed before i.e. SN la, BAO, CMB and WL, a number of other 
methods have been proposed, giving the possibility to cross check the same parameters. We 
report here the most promising ones: 

• Gamma ray bursts: GRB are powerful explosions, probably generated by the collapse 
of a rapidly-rotating, high-mass star into a black hole. They have a huge luminosity. 
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Figure 17: Constrains on {ilm,crs) from weak lensing observations CFHTLS and the CMB 
observation WMAP3 adopting ACDM. 




Figure 18: Confidence region at 68% for the dark energy parameters for SNAP experiment, 
which combines 2000 SNe to z = 1.7 and a weak lensing survey of 1000 sq. deg. Left panel: 
Constraints in the ^l,n-w plane, assuming constant w. Right panel: Constraints in the wq- 
Wa plane for time-varying dark energy equation of state, marginalized over ri,„ for a flat 
Universe. 
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which can be seen at very large distances. They could become the competitor of SN 
la as they can cover different redshifts (up to z = 8.3, see Ref. |257) ) hence testing the 
behavior of the dark energy in different epochs. Although the physics of the GRB is 
still uncertain, several papers have been published in order to find a correlation between 
the luminosity of GRB and other observables so that they could be used as standard 
candles. So far the constraints on dark energy parameter are very weak. 

• Integrated Sachs- Wolfe effect: when photon enters a potential well it gains energy while 
climbing out of a potential loses energy. During matter domination era the gravitational 
potentials are constant, hence a photon falling down and climbing up will not manifest 
a change in the energy. If matter is not the dominant component, the gravitational 
potentials are expected to evolve. As a consequence, photons traveling through a 
shallowing (or deepening) potential well will manifest a blueshift (or redshift). This is 
known as ISW effect. We have shown in the previous sections that the dark energy 
affects the evolution of the gravitational potentials either at background level and at 
the perturbative level. At the present the ISW effect has been measured up to redshift 
z « 1.5, see Refs [25811259] . 

• Age tests: the Hubble parameter is defined as H = —j^j^', if we know the age 
difference At of two galaxies separated by the redshift Az then it is possible to measure 
the Hubble parameter H{z) directly, see Ref. IgSOl UMl IISI UMl HSi] . The idea 
of using stellar age seems very interesting since it will use a complete different method. 
Here we are not interested in geometry but on cosmic chronology. Of course the difficult 
part of this method is that we have to make sure that the clocks we are using are reliable. 

• Gas fraction: in standard cosmology is expected that all clusters contain the fixed ratio 
of baryons to total matter ^b/^m = const. Now, the X-ray thermal bremsstrahlung 
luminosity that comes from those baryons present in galaxies and clusters of galaxies is 
proportional to the volume of the emitting region which is proportional to the total gas 
mass. It follows that the gas mass is proportional to the luminosity. An observer will 
measure the fiux which is a function of the luminosity and of the luminosity distance 
which depends on the cosmology. The constraints obtained using these methods, see 
Refs. [2661 12671 1268j . are still not accurate because of the different assumptions made, 
like spherical model, ideal gas equation of state, isothermal distribution, etc. 

12 Conclusion 

There is now strong evidence that the Universe is accelerating. However the main question is: 
is it really accelerating or it is just an observed effect probably due to other unknown physics? 
And, if it is really accelerating, is this acceleration due to a strange, small component called 
dark energy? Finding out the real nature of the dark energy is one of the most exciting and 
challenging problems facing physicists. In this review we have not investigated all the dark 
energy models because there have been more than 1000 papers with the dark energy in the 
title written since 1998. In this review we were more interested in the possible explanation 
of the acceleration of the Universe focusing of course on dark energy but also highlighting 
alternative solutions. 

After more than 10 years, dark energy still remains an open question. What we have 
done is firstly to add a new fluid component in the Universe, concentrating on the possibility 
that the dark energy may be dynamical; furthermore we also allowed the possibility that 
the Einstein equations themselves require modification. Although there is not yet reason 
to do so, we think it was appropriate to consider alternative theories which in principle are 
plausible. 



62 



We paid special attention to the observational tests because the list of planned and 
imdergoing observational projects related to dark energy is impressive, see Ref. |269) . The 
amount of data we will possess in the future will be enormous and we will be able to constraint 
other parameters apart from the classical cosmological ones. This will help us to say a bit 
more about the candidate responsible for the expansion of the Universe. It might be that 
we need to go beyond the cosmological constant A or it is also possible that the experiments 
will keep confirming A as the best candidate. In the last case, we will learn more about 
gravity or we will be able to map the total matter content or we will be able to justify all the 
assumptions we have been making so far. For instance, it has been shown via weak lensing 
and galaxy velocity that General Relativity is confirmed also at large scales, see Ref. |270) . 

Again, this review should not be seen as a list of all the models on dark energy and the 
possible suggestions to explain several problems. At the moment, it is impossible to say 
anything! In other words, the term Dark Energy is the name that we give to something that 
we still do not understand to explain the observed late time acceleration of the Universe. 

Moreover, this review was not meant to discourage the reader because she (or he) might 
have thought that a lot was already done in order to understand the Universe where we live. 
On the contrary, this review shows that there is still a lot to do. 

In the end, studying and understanding the dark energy problem will not be a waste of 
time and money but it will most likely be an access to a complete new world. 
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